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Closed  form  solutions  to  the  stress  analysis  for  an  elastic  coil  of  superconducting 
magnets,  which  include  isotropic  Green's  functions  solution  and  orthotropic  generalized 
plane  strain  analysis,  are  presented.  The  prediction  of  stress  and  strain  is  essential  for  both 
mechanical  and  electrical  design  of  high  field  superconducting  magnets.  The  generalized 
plane  strain  analysis  provides  an  analytical  formulation  for  the  axisymmetric  state  of  stress 
and  strain  in  the  midplane  of  a  superconducting  magnet,  where  shear  stress  is  negligible. 
This  formulation  takes  into  account  the  effect  of  the  axial  body  force.  The  full 
axisymmetric  stress  analysis  in  terms  of  Green's  functions  provides  the  components  of 
stress  throughout  the  coil  and  includes  the  shear  stress  in  addition  to  the  normal  stresses. 
Green's  functions  method  permits  the  development  of  a  solution  irrespective  of  the  type  of 
the  field  or  its  distribution  within  the  coil.  Green's  functions  are  derived  by  using  finite 
integral  transforms  in  an  interval  appropriate  for  a  coil,  e.g.,  finite  Hankel  transform.  A 
comparison  of  the  analytical  results  with  the  finite  element  methods  is  presented  to  confirm 
the  validity  of  the  solutions.  The  results  show  good  agreement  for  both  generalized  plane 
strain  analysis  and  Green's  functions  solution  with  the  results  of  finite  element  methods. 
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CHAPTER  1 
INTRODUCTION 


Stress  analysis  for  a  high  field  solenoid  magnet  is  an  important  aspect  of  the  overall 
magnet  design.  The  prediction  of  stress  and  strain  is  essential  for  both  mechanical  and 
electrical  design  of  high  field  solenoid  magnets.  These  magnets  are  designed  in  a  variety  of 
configurations.  A  superconducting  magnet  is  one  example  of  such  magnets,  which  can  be 
treated  as  a  combination  of  several  coils,  where  each  coil  may  be  reinforced  by  a 
nonconducting  layer.  Depending  upon  the  geometrical  specifications  of  a  coil,  magnetic 
fields  may  behave  differently.  These  fields  result  in  magnetic  body  forces,  and  thus 
stresses,  in  all  directions.  Traditionally  only  the  tangential  component  of  stress  has  been 
considered  for  design  and  failure  analysis.  Recent  experimental  results  indicate  that  in 
addition  to  the  normal  stresses,  shear  stress  also  contributes  to  failure.  The  value  of  shear 
stress  has  been  determined  to  be  small  in  the  midplane  and  becomes  larger  toward  the  ends 
of  the  coil.  In  the  present  analytical  solutions,  shear  stress  is  assumed  to  be  negligible  and 
the  use  of  finite  element  methods  is  time  consuming  especially  for  the  case  of  designing  a 
magnet  with  many  layers.  As  a  result  a  three-dimensional  closed  form  solution  is  desired 
to  develop  design  criteria  and  understand  the  distribution  of  stresses  (including  shear) 
through  a  solenoid  magnet.  Such  a  closed  form  solution  can  be  used  in  a  numerical  scheme 
for  multilayered  magnet  structures. 

In  this  work,  the  general  closed  form  solutions,  using  Green's  function  method  and 
generalized  plane  strain  condition,  are  derived  for  high  field  solenoid  magnets.  These 
solutions  are  applied  to  the  important  case  of  superconducting  magnets.  The  generalized 
plane  strain  analysis  is  performed  for  the  midplane  of  an  elastic  orthotropic  coil  of  a 
solenoid  magnet  where  shear  stress  is  negligible.  The  full  axisymmetric  stress  analysis  for 
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an  elastic  isotropic  coil  is  derived  using  the  Green's  functions  solution.  This  analysis  is  not 
limited  to  the  midplane  and  can  be  used  for  any  type  of  solenoid  magnet. 

Historically,  magnet  stress  analysis  has  been  presented  as  a  plane  strain  or  plane 
stress  problem.  Here,  a  better  treatment  in  terms  of  generalized  plane  strain  is  developed, 
which  takes  into  account  the  effect  of  the  axial  body  force.  The  generalized  plane  strain 
provides  an  accurate  prediction  of  all  normal  components  of  stress  and  as  a  result  is  a  good 
approximation  to  the  general  three-dimensional  solution.  For  the  first  time  a  three- 
dimensional  stress  analysis  for  an  isotropic  axisymmetric  coil  with  finite  intervals  in  both 
radial  and  axial  directions  and  subjected  to  body  forces  is  derived.  The  use  of  Green's 
functions  is  the  best  possible  technique  for  providing  solutions  for  any  type  of  body  forces. 
In  this  analysis,  the  boundary  conditions  for  the  coil  become  complicated  and  cannot  be 
satisfied  by  the  traditional  methods.  A  new  approach  is  used  to  satisfy  the  boundary 
conditions  which  results  in  a  complementary  solution. 

The  use  of  Green's  functions  solution  is  not  limited  to  magnetic  body  forces.  It  can 
be  applied  to  other  axisymmetric  elasticity  problems  with  finite  intervals.  The  Green's 
functions  solution  can  be  used  for  inclusion  problem  in  composite  materials  where 
eigenstrains  may  be  considered  as  body  forces.1"4  The  Green's  functions  solution  can  also 
be  applied  to  crack  problems,  due  to  residual  stresses,  in  composite  materials,  where 
difference  between  the  thermal  expansion  coefficients  of  the  fiber  and  the  matrix  results  in 
fictitious  body  forces.5-6 

This  dissertation  is  divided  into  six  chapters  and  five  appendices.  In  Chapter  two 
superconductivity  is  briefly  explained  and  magnetic  fields  for  a  circular  wire  are  derived. 
In  Chapter  three  magnetic  and  thermal  stress  analyses  based  upon  generalized  plane  strain 
assumption  for  an  elastic  orthotropic  coil  are  introduced  and  the  arbitrary  constants  in  the 
magnetic  stress  formulation  are  determined  for  different  kinds  of  boundary  conditions.  The 
derivation  of  the  arbitrary  constants  for  the  thermal  stress  formulation  is  included  in 
Appendix  A. 
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The  main  theme  of  this  work  is  in  Chapter  four,  where  the  three-dimensional 
axisymmetric  stress  analysis  in  terms  of  Green's  functions  is  performed.  The  problem  is 
formulated  in  terms  of  radial  and  axial  stress  functions.  The  partial  differential  equations 
for  stress  functions  are  derived  from  the  fundamental  field  equations  (a  detailed  derivation 
of  the  partial  differential  equations  is  given  in  Appendix  B).  The  radial  stress  function  is 
solved  in  terms  of  the  radial  Green's  function  by  using  the  finite  Hankel  transform  of  order 
one  and  the  finite  Fourier  cosine  transform.  The  axial  stress  function  is  derived  in  terms  of 
the  axial  Green's  function  by  using  the  finite  Hankel  transform  of  order  zero  and  finite 
Fourier  sine  transform.  The  derivation  of  the  finite  Hankel  transform  is  also  included  in 
this  chapter  (proof  of  orthogonality  property  of  the  finite  Hankel  transform  is  provided  in 
Appendix  C).  The  boundary  conditions  are  satisfied  by  a  complementary  solution  for  the 
axial  stress  function.  By  using  a  new  approach,  the  complementary  solution  is  derived  for 
traction  free  boundary  conditions.  The  derivation  of  the  complementary  solution  is  not 
limited  to  the  traction  free  boundary  conditions  and  can  be  determined  for  any  kinds  of 
boundary  conditions.  The  final  solution  for  stresses  is  obtained  from  the  radial  Green's 
function  together  with  the  superposition  of  the  axial  Green's  function  and  the 
complementary  solution.  Moreover  in  this  chapter,  for  each  of  the  derivations  of  similar 
equations  only  one  derivation  of  an  equation  is  given  and  the  rest  are  provided  in 
appendices  (Appendix  D,  Appendix  E  and  Appendix  F). 

In  Chapter  five  the  results  of  stress  analysis  for  the  example  magnets  are  presented. 
These  results  are  compared  to  the  finite  element  methods  to  confirm  the  validity  of  the 
analytical  stress  analysis.  In  Chapter  six  the  conclusion  is  provided  and  some  future  work 
is  proposed. 


CHAPTER  2 
SUPERCONDUCTING  MAGNETS 


Introduction 

A  superconducting  magnet  is  a  nested  set  of  individual  cylindrical  superconducting 
coils  on  separate  support  structures  (Fig.  2.1).  A  superconducting  coil  may  have  additional 
support  structure  (reinforcement)  on  the  outer  diameter  to  constrain  the  coil.  There  is  a 


Figure  2. 1  Schematic  diagram  of  a  representative  superconducting  magnet. 

long  history  of  attempts  to  understand  and  improve  the  performance  of  superconducting 
magnets.  The  tendency  of  coils  to  quench  prematurely  at  relatively  low  fractions  of  the 
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critical  current  is  often  attributed  to  mechanical  issues.  Stress  analysis,  and  a  knowledge  of 
stress,  strain  and  displacement  of  the  windings,  is  therefore  essential  to  the  design  of 
superconducting  magnets. 

Superconductivity 

The  complete  absence  of  electrical  resistivity  for  the  passage  of  current  below  a 
certain  "critical"  temperature  (less  than  21  °K)  is  the  basic  premise  of  superconductivity.  In 
addition  to  critical  temperature,  the  critical  field  and  critical  current  density  are  two 
parameters  that  define  a  critical  surface.  Fig.  2.2  shows  a  surface  in  a  three  dimensional 
space  whose  axes  are  temperature  T,  magnetic  field  B  and  current  density  J. 

I  J 


B 


T 


Figure  2.2  Critical  surface  of  a  typical  superconducting  magnet. 

The  superconducting  phase  exists  only  for  the  points  below  the  surface;  points  above  the 
surface  correspond  to  the  presence  of  resistivity.8  The  critical  temperature,  Tc,  of  a  magnet 
is  the  temperature  where  resistivity  suddenly  drops  from  a  high  value  (10  uXXcm)  to  a  low 
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value  (less  than  10"9  u\Q.cm).  For  comparison,  resistivity  of  the  purest  copper  is  about 
10"3  p-Q.cm  so  for  all  practical  applications  resistance  of  a  superconductor  is  assumed  to  be 
zero.  The  transition  between  normal  state  and  superconducting  phase  takes  place  at  the 
critical  field,  Bc.9  The  critical  current  density,  Jc,  is  the  maximum  current  density  that  can 
be  applied  to  a  conductor  in  the  superconducting  phase.  The  critical  current  density  of  a 
superconducting  wire  is  strain  dependent.10  This  is  specially  true  for  Nb3Sn 
superconductors  as  shown  in  Fig.  2.3.  For  most  superconducting  wires,  the  critical 
current  is  observed  to  increase  initially  with  tensile  strain  and  then  decrease  with  increasing 
tension. 


T3 

N 


o 
Z 


0.6  0.8 
Strain  % 
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1.2 


Figure  2.3  Uniaxial  strain  dependence  of  the  critical  current  density  of  a  Nb3Sn 
superconductor.8 


Magnetic  Fields 

The  state  of  stress  in  a  coil  subjected  to  a  uniform  current  density  is  a  direct  result  of 
magnetic  fields.  Magnetic  fields  for  a  solenoid  coil  may  be  calculated  by  integrating  the 
contributions  from  individual  circular  current  filaments.  The  magnetic  field  vector,  B, 
from  a  circular  current  loop  is  calculated  from  the  vector  potential  A,  where  B  =  V  x  A . 
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The  vector  potential  A  at  point  P  (see  Fig.  2.4)  due  to  a  current  /  has  only  a  tangential 
component  given  by  Eq.  (2.1). 11 

where  /x0  is  the  magnetic  permeability,  a  is  the  radius  of  the  loop  and  R  is  a  vector  from 
point  P  to  any  point  on  the  loop.  Vector  R  may  be  written  in  the  form  R  =  oP  -  oB.  In 
cylindrical  coordinates  (r,  (p,  z),  oB  and  oP  are  expressed  by  acos<joer  +  asin<pe  and 
rer  +  zez,  respectively.  Using  these  expressions,  we  may  write 

R  =  (r-acos<p)er  -asin^^  +  zez 

and,  hence, 

|R|  =  ^(r-acos(p)2  +a2sin2  cp  +  z2 

or  by  expanding  the  (r  -  a  cos  (p)2 , 

|R|  =  ^r2  +a2  +z2  -2arcos<p .  (2.2) 

The  closed  path  of  the  integral  in  Eq.  (2.1)  is  a  circle,  thus,  ds  =  ad(p.  Introducing 
Eq.  (2.2)  into  Eq.  (2.1)  gives 


_  fi0I  r  acosq>d(p 

\(r'Z)-  ~T~  J      I    7  1  9 

47r  -VVr  +a  +z  -  2arcos (p 


or  since  the  expression  inside  the  integral  is  an  even  function  in  terms  of  <p 

A<       M0/f  a  cos  (pd(p 

AP<r'rt-—    I   /  ,  9 — ,  (2.3) 
27r  o  Vr  +a  +z  ~2arcos(p 


The  substitution  of  (p  =  n-  29  into  Eq.  (2.3)  yields  Eq.  (2.4). 
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acos29d6 


n  {  ^r2+a2+z :  +  2arcos20 


(2.4) 


Figure  2.4  Magnetic  field  at  any  point  in  space  from  an  ideal  circular  loop. 


By  substituting  cos 20  =  1  -  2 sin2  6  in  Eq.  (2.4),  we  may  write 


Kj2 


fU  \  a(l-2sin20)J0 


x  {  V(r  +  a)2  +Z2  -4ar  sin2  9 


(2.5) 


Introducing  a  new  parameter  it  = 


2^[ar 


V(r  +  a)2+z2 


into  Eq.  (2.5),  yields 


__VoI  1 


71  4{r  +  a)2+z' 


fad-: 


a(l  -  2 sin2  fl)Jfl 
it2  sin2  6  ' 
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or 


Am(r,z)  =  - 


V         a  7  dO 

Tt  ^(r  +  af+z2  [I  Vl-*2sin20 

2 


(2.6) 


Y(l-fe2sin20)j0  * 
{   Vl-fc2sin20  l 


fc2sin20 


The  integrals  in  Eq.  (2.6)  are  in  form  of  elliptic  integrals.  Thus,  we  may  write 


V-*>  =  ~—  ;      a2     ,  {g«)  +  ^-[Ew-gW]}  (2.7) 


where  A'ajand  are  elliptic  integrals  of  the  first  and  second  kind  given  by  Eqs.  (2.8) 
and  (2.9),  respectively. 

*/2 


A"(*)  =  | 


de 


Vl-jt2  sin2  0 


(2.8) 


*/2 

£(*>=  jVl  -Jfc2sin2 


(2.9) 


The  substitution  of 


^(r  +  af+z2  lja~r 


into  Eq.  (2.9),  yields 


'  fc2^ 
1-  — 

V  2y 


K(k)-E(k)\ 


(2.10) 


The  components  of  magnetic  field  are  calculated  from  B  =  VxA,  where 


A  =  A„e9. 


Or(r. z)  —  -  7=  ====■<  —A (t 

A      2k  r  V(r  +  fl)2  +z2  I 


2.2        2  1 

<-  r  +a2  +  Z2  E(k)\  (2.11) 


(r-a)2+z 
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=  -^-(rAf)  =  M         1         L(t)  +fl2    ^/agJ  (2.12) 
r<?r  2;r  ^r  +  af  +  z2  {         (r-a)2+z2  J 

Eqs  (2.11)  and  (2.12)  indicate  that  the  axial  and  radial  fields  are  even  and  odd 
functions  of  z,  respectively 


CHAPTER  3 
GENERALIZED  PLANE  STRESS  ANALYSIS 

Introduction 

An  axisymmetric  stress  analysis  of  superconducting  magnets  is  formulated  under 
assumptions  that  make  it  applicable  to  a  region  about  the  midplane  of  a  long  magnet.  The 
stress  analysis  of  superconducting  magnets,  as  an  important  aspect  of  the  overall  design 
process,  has  been  the  subject  of  a  number  of  publications.12-19  Early  analysis  assumed 
only  isotropic  material  properties  and  plane  conditions,  and  treated  only  the  dominant 
tangential  stress.12"14  These  analyses  were  further  developed  to  include  orthotropy  and 
solutions  to  the  winding  and  thermal  stresses.15-19 

The  generalized  plane  strain  analysis  provides  an  analytical  formulation  of  the 
axisymmetric  state  of  stress  and  strain  in  the  midplane  of  a  superconducting  magnet.  The 
present  formulation  takes  into  account  the  effect  of  the  axial  body  force  between  the 
midplane  and  the  end  of  the  coil.  A  fundamental  simplifying  assumption  is  that  the  shear 
stress  in  the  midplane  is  zero.  A  second  assumption  which  leads  to  the  present  solution  is 
that  the  axial  strain  in  the  midplane  region  of  a  long  magnet  is  constant  through  the  coil 
along  a  radius,  for  a  fixed  axial  location.  The  applicability  of  this  assumption  for  the  stress 
analysis  of  superconducting  magnets  has  been  recognized  previously.20- 21 

Magnetic  Stress  Analysis 

Stress  and  Strain  Formulation 

The  magnetic  stress  and  associated  strain  are  the  results  of  the  distributed  Lorentz 
forces.  The  equations  for  the  magnetic  stress  and  strain  are  formulated  for  a  coil  of  a 
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12 


superconducting  magnet  with  a  uniform  current  density.  The  coils  of  a  superconducting 
magnet  are  treated  as  homogeneous,  orthotropic  and  linear  elastic  material. 

With  the  assumption  of  zero  shear  stress,  the  constitutive  equations  relating  stress 
and  strain  are  given  by  Eq.  (3.1). 

£e=te~Vret~Vz8t  (3'1} 
Ee    Er      w  Ez 

Oft  <Xr  C\ 

£=-v*r  —  -vr7  —  +  —L 
ZEe      nEr  Ez 

Eq.  (3.1)  may  be  written  in  the  form  of  the  associated  compliance  matrix  S. 

e,=S,joj  (3.2) 


Inverting  Eq.  (3.2),  the  stress  and  strain  are  related  by  the  stiffness  matrix  C. 

d  =  dj£j 


(3.3) 


Following  the  traditional  approach,  the  radial  displacement  u  is  introduced  as  a 
primary  variable.  The  tangential  and  radial  strains  are  related  to  the  displacement  as  given 
in  Eq.  (3.4). 

u  du  A^ 

£r  =  -r  (3.4) 
r  dr 

In  Eq.  (3.5),  the  components  of  the  stress  are  written  in  terms  of  the  displacement  and  the 
axial  strain. 

°e  =  Cu-  +  Cn  —  +  C,3£J  (3-5) 
r  dr 

u     _  du 

Or  =  C\2~+  Cn—  +  C2i£z 

r  dr 

u    ^  du 
<Jz  =  Ci  3  -  +  C23  —  +  C33  £z 
r  dr 
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It  is  assumed  that  the  axial  strain  is  constant  throughout  the  coil. 

£,  =  constant  (3.6) 

In  the  absence  of  the  shear  stresses,  the  equilibrium  equations  at  the  midplane  are 
reduced  to  one  equation  as  given  by  Eq.  (3.7). 

j-{ror)-<Je  +  rXr=0  (3.7) 

The  radial  body  force  Xr  is  due  to  the  radial  Lorentz  force  density  given  by 

Xr  =  /,flt(r)  (3.8) 

where  Je  is  the  current  density  and  2?,(o  is  the  axial  magnetic  field. 
Introducing  Eqs.  (3.5)  and  (3.8)  into  (3.7)  gives 

(  du\ 


Q2  , 
dr 


r  — 

\  dr ) 


~  Cu  " +(Q3  -  Cl3)ez  =  -rJeBzir).  (3.9) 


It  is  assumed  that  the  radial  distribution  of  the  axial  field  is  linear. 

Bzir)  =  Bc-C0r  (3.10) 

The  constants  Bc  and  C0  are  determined  from  the  values  of  the  magnetic  fields  at  the  inside 
and  outside  radii.  The  accuracy  of  the  linearity  assumption  has  been  previously  studied  in 
some  detail.18  This  assumption  is  not  fundamental  to  the  method,  and  with  a  few  changes 
in  the  algebra  a  high  order  polynomial  for  the  field  may  be  adopted. 

Incorporating  Eq.  (3.10)  in  Eq.  (3.9)  yields  Eq.  (3.11),  where  the  customary 
variable  k  is  the  anisotropy  factor  given  in  Eq.  (3.12). 

d2u    du    ,2u       J0Br      JRCn  2  (C„-C.,) 
r — -  +  k2-  =  +        r2  (3.H) 

dr2     dr        r        C22        C22  C22        z  K  ) 
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The  general  solution  of  Eq.  (3. 1 1)  may  be  written  in  the  form 

u  =  Dxrk  +D2r~k  +  Aiezr  +  A2r2  +A,r3  (3  13) 

where  D\  and  are  integration  constants,  and  A,,  A1  and  A3  are  constants  based  upon 
the  coil's  magnetic  and  mechanical  properties. 

A=-(i-*2)    c22  (3-14) 

(4-k2)  C22 
1      /  C 


4  = 


(9-*2)  C 


22 


By  inspection,  there  are  three  singular  values  of  the  anisotropy  factor  k.  When  k  has  the 
value  1,  2,  or  3,  the  associated  constant  A  is  replaced  in  the  general  solution  by  the 
corresponding  function  A  given  by  Eq.  (3.15),  (3.16),  or  (3.17). 


r\  _      1  (Q3 


A,(r)  =  --^  ^lnr  (3.15) 


2  C 


A2(r)  =  -i^lnr  (3.16) 
4  C22 


A3(r)  =  +I^.lnr  (3  17) 


6  C22 


In  the  remainder  of  the  analysis,  non-singular  values  of  k  are  assumed. 

The  solution  for  the  displacement  is  used  in  Eqs.  (3.4)  and  (3.5)  to  yield 
expressions  for  the  state  of  strain  and  stress  as  given  in  Eqs.  (3.18)  through  (3.22), 
together  with  Eq.  (3.6). 


£e  =  D/_I  +  D2r~k-'  +Alez+A2r  +  A3r2 


(3.18) 
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£T  =  £D,r*~'  -  kD2r'k'x  +Alez+  2A2r  +  3A3r2  (3.19) 

ae  =  (C„  +kCn)D/-1  +(C„  -kCn)D2r-k-x  (3.20) 
+(C„  +Ci2)Ai£z  +(C„  +  2Cl2)A2r  +  (Cu+3Cn)A3r2  +  Cl3£z 

ar  =  (Cu+kC22)D/-'+(Cl2-kC22)D2r-k-'  (3.21) 
+(C12  +  C22)A£,  +  (C12  +  2C22)  V  +  (Cl2  +  3C22)  V  +  C23£2 

crz  =  {Cl3  +  kC23)D/-l+{Cn-kC23)D2r-k-1  (3.22) 
+(C|3  +  C23 

)V«  +  (C13  +  2C23)V  +  (C13  +  3C23  )A3r2  + 

These  equations,  Eqs.  (3.18)  -  (3.22),  together  with  the  values  of  the  integration  constants 
D  i  and  D2  plus  the  value  of  the  axial  strain  £z ,  determine  the  distribution  of  stress  and 
strain  in  a  given  coil. 

Equations  for  the  Constants 

Solutions  for  the  integration  constants  D  \  and  D2  results  from  the  boundary 
conditions  associated  with  each  radial  layer.  For  a  stand-alone  coil  with  a  single  layer,  the 
radial  stress  at  the  inside  and  outside  radii  will  be  zero.  For  a  coil  with  several  radial 
layers,  additional  conditions  are  extracted  from  the  continuity  of  the  radial  stress  and  radial 
displacement  at  the  interface  between  each  layer.  In  this  way,  the  number  of  equations 
match  the  number  of  integration  constants.  An  additional  equation  is  required  to  determine 
the  unknown  value  of  the  axial  strain. 

The  basic  assumption  is  made  that  at  a  given  axial  position,  e,  is  constant  in  each 
layer,  and  that  for  coils  with  more  than  one  layer,  the  value  of  the  axial  strain  is  the  same  in 
all  layers,  including  external  reinforcement.  This  assumption  is  equivalent  to  the 
assumption  that  the  coil  remains  a  right  cylinder  in  a  neighborhood  of  the  axial  plane  in 
which  the  stress  is  calculated. 
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Using  the  concept  of  a  plane  through  the  coil  at  a  given  axial  location,  the  static 
equilibrium  of  the  coil  requires  that  the  net  axial  stress  over  the  plane  is  equal  to  the  total 
applied  axial  load  accumulated  on  the  plane.  Thus,  on  a  plane  through  the  coil  at  an  axial 
position  z, 

\lnrazdr  =  Fz  (3.23) 

where  a\  and  an+\  are  the  inside  and  outside  radii  of  a  coil  which  has  n  distinct  radial 
layers,  and  Fz  is  the  total  axial  Lorentz  force  between  the  midplane  and  the  end  of  the  coil 
over  all  layers. 

For  a  single  layer  coil  the  boundary  conditions  are 

<Tr=0  at  r  =  fl,  (3.24) 
CTr  =  0  at  r  =  a2  (3.25) 

\l7trOzdr  =  Fz. 

Evaluating  Eq.  (3.21)  at  the  inner  radius  (ai)  results  in 

auD,+anD2+al3£z=bl 

where  the  constants  are  given  by 

«„  =(C12  +  fcC22)af-' 
0,2  =  (C12  —  kC22)ax 

°13  =  (^12  +  +  Q3 

h  =  "(C12  +  2C22)A2a1  -  (da  +  3C22)V,2- 
Evaluating  Eq.  (3.21)  at  the  outer  radius  (ai)  leads  to 


(3.26) 


(3.27) 


(3.28) 
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(3.29) 


where  the  constants  are  given  by 

a22  =  \C\2  ~  kC22)a 

°23  =  aU 

l?2  — ~~*\^12      ^^22  ),^2^2  ~~ '  (^12      ^^22  )^3^' 


-k-\ 


(3.30) 


Integrating  Eq.  (3.22)  in  Eq.  (3.26)  gives 


a^Dx+anD2+a^ez=b^ 


(3.31) 


where  the  constants  are  given  by 

a3i  =  (^13  +  ^23) 
a32  =  (Ci3  — 


k  +  l 

-k+\  -k+l 
a2  ~a\ 


-k  +  l 

«33=[(C13  +  Q3)A+C33] 


„2  2 

a2  —  a, 


(3.32) 


«2  - 


The  set  of  linear  Eqs.  (3.27),  (3.29),  and  (3.31)  may  be  written  in  matrix  form 


an 

"l3 

A 

V 

<h\ 

a22 

"23 

A 

fl3l 

«32 

"33. 

(3.33) 


which  is  solved  in  the  standard  way  as  given  by  Eq.  (3.34). 


au 

an 

«13_ 

-1 

y 

A 

a2\ 

a22 

«23 

b2 

.«31 

"32 

"33. 

A. 

(3.34) 
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For  a  single  layer  coil  with  a  reinforcement,  the  boundary  conditions  are  applied  to 
the  coil  and  the  reinforcement. 


=  0  at  r  =  ax  (3.35) 


G(rl)  =  a(r2)  at  r  =  a2  (3.36) 


M<°  =  w<2)  at  r  =  a2  (3.37) 


cr<2)  =0atr  =  a3 


8| 


^2Kra.dr  =  Fz 


Evaluating  Eq.  (3.21)  at  the  inner  radius  (a\)  results  in 


auDl+al2D2+al5ez=bx 


(3.38) 


(3.39) 


(3.40) 


where  the  constants  are  given  by 

an={Cn+kC12)al'  (3.41) 
an  ~  (C12  -  kC22)ax 
a,5  =  (C12  +  C22)A,  +  C23 
^.  =-{Cn+2C22)A2al-(Cl2+3C22)A3al 


Evaluating  Eq.  (3.21)  at  the  interface  between  the  coil  and  the  reinforcement  (ai)  results  in 

a2x Dx  +  a22D2  +  a2iD[  +  a24D2  +  a25sz  =b2  (3  42) 

where  the  constants  are  given  by 
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<hi  ={Cn+kC22)ak2-1  (3.43) 

a2i  =  ~ \C\2  ~*~  ^  ^2l)a2 

a24  =  —[Cn—k  C22)a2 

a25  =  [(^12  +  Q2  ) A  +  Q»  ]  ~  [(^12  +  Q'2  )  A'  +  ^23  ] 

Z?2  =  — (C12  +  2C22)A2a2  —  (C12  +  3C22)A3a2 


-t'-i 

*2 


and  the  unprimed  and  primed  quantities  refer  to  the  coil  and  the  reinforcement,  respectively. 
Evaluating  Eq.  (3.13)  for  the  displacement  at  the  interface  results  in 

<hA  +  fl32  D2  +  fl33Dl'  +  «34  A'  +  ai5£z  =  bl  (3.44) 

where  the  constants  are  given  by 

«3,=fl2*  (3.45) 


a32  =  a2 


«33  =  "«2 


t' 
2 

-*' 


«34  =  -«2 
a35  =  Afl2  -  Afl2 

Z>3  =  - A,a2  -  A3a2 


Evaluating  Eq.  (3.21)  at  the  outside  radius  of  the  reinforcement  (03)  results  in 

a4,D;  +  a44D^  +  a45ez=0 


(3.46) 


where  the  constants  are  given  by 

fl43  =  (^12  ^-"22)^3  ^  ^ 

a44  =  (^12  —  k  ^"22)^3 

a45  =  (^12  ~*~  ^22)  A  +  ^23' 

Integrating  Eq.  (3.22)  through  coil  and  reinforcement  in  Eq.  (3.39)  results  in 
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"51  A  +  "52  A  +  "53  A'  +  "54  A'  +  "55CZ  =  ^5 


(3.48) 


where  the  constants  are  given  by 
"51  =  (^-13  +  ^23) 


k  +  \  k+l 
a2  ~"l 

k  +  l 

„-k+\ 

a2     —  a, 
-k  +  l 


jt'+i  _  t'+i 


"52  —  (^13      ^23  J 


k'  +  l 


k'+\  -k'+\ 


"s4  ~  (^13      ^  Q3) 


\  "3  ~a 


-k'  +  \ 


"55  =  [(c13  +  c23  )a,  +  c33  ]       +[{c;3  +  c2'3  )a; + q3  ] 


2  _  2 

\<h-<h 


3  3 

a,  —  a, 


b5  =  -(C13  +  2C23)^  _  (Ci3  +  3C32  )A3  +  £. 


(3.49) 


The  set  of  linear  Eqs.  (3.40),  (3.42),  (3.44),  (3.46),  and  (3.48)  may  be  written  in  matrix 
form 


"u 

a\2 

0 

0 

"is" 

"A" 

V 

a21 

a22 

"23 

a24 

A 

*2 

"31 

«32 

fl33 

a34 

fl35 

A 

0 

0 

«43 

"44 

°45 

A 

h 

"51 

fl52 

"53 

a54 

fl55. 

A  j 

A. 

(3.50) 


which  is  solved  as  expressed  by  Eq.  (3.51). 


"A" 

a12 

0 

0 

"l5_ 

-1 

V 

A 

a22 

a23 

"24 

"25 

^2 

A 

"3I 

"32 

"33 

"34 

"35 

*3 

A 

0 

0 

"43 

"44 

"45 

^4 

A. 

a5I 

"52 

"53 

"54 

a55. 

A. 

(3.51) 


For  a  general  multilayer  coil  («  layers)  the  boundary  conditions  are  applied  to  each 


layer. 
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cx'0  =0  at  r  =  a. 


CT(;)  =  (T(;+i) 


«ry)  =  ««+1)  at  r  =  ay+I;j  =  U-l 


at  r  =  aj+l;  j  =  l,n-l 


<J™  =  0  at  r  =  an+i 


(3.52) 
(3.53) 
(3.54) 

(3.55) 


The  axial  force  equilibrium  also  applies  to  all  layers. 

j2Kr(Jzdr  =  Fz 


(3.56) 


The  derivation  of  the  equations  for  the  constants  proceeds  as  in  the  previous  cases. 
The  coefficient  matrix  for  the  system  of  2n+l  equations  and  unknowns  is  given  by  Eq. 
(3.57). 


an 

0 

0 

0 

0 

0 

an 

<hi 

0 

0 

0 

<hi 

an 

a33 

a34 

0 

0 

0 

0 

0 

"43 

«46 

0 

0 

0 

"54 

"56 

0 

0 

0 

1  2«  +  l 


2  2n+\ 


a 


3  2«+l 


*4  2n  +  l 


a 


5  2n+l 


^       fl2/i  2n-l      a2«  2n 


2«  2n  +  l 


(3.57) 


_a2«+I  I      a2«  +  12      '*'  a2n+\2n      'hn  +  l  2n  +  l 

The  solution  for  this  system  of  equations  is  obtained  by  matrix  inversion. 

Thermal  Stresses  Analysis 

The  thermal  stress  and  associated  strain  are  the  result  of  anisotropy  in  thermal 
contraction  coefficients  within  a  coil,  or  differences  in  thermal  contraction  coefficients 
between  the  layers  including  reinforcement.  The  equations  for  thermal  stress  and  strain  are 
formulated  for  a  coil  of  a  superconducting  magnet. 


22 


With  the  assumption  of  zero  shear,  the  stress  is  related  to  the  mechanical  component 
of  the  total  strain  by  the  stiffness  matrix  C. 

<Ji  =  C'J£j  (3.58) 

It  is  recognized,  in  dealing  with  change  in  temperature,  that  the  total  strain  within  a  body  is 
the  sum  of  strains  produced  by  mechanical  stresses  and  strains  associated  with  thermal 
expansion  or  contraction.  The  radial  displacement  u  ,  that  is  introduced  as  a  primary 
variable  in  the  analysis,  is  directly  related  to  the  total  strains  in  the  tangential  and  radial 
directions.  The  mechanical  strains  are  given  by 

(3.59) 


_  u 

■  CCgAT 

r 

du 

-arAT 

e, 

~~dr 

-azAT 

where  e'°'  is  the  total  axial  strain.  Introducing  Eq.(3.59)  into  Eq.  (3.58)  gives  the  stress  in 
terms  of  the  displacement  and  the  total  axial  strain.  The  generalized  plane  strain  assumption 
applies  to  the  total  axial  strain. 


£'"'  =  constant  (3.60) 


In  the  absence  of  the  shear  stresses  and  body  forces,  the  equilibrium  equations  are 
reduced  to  one  equation  as  defined  by  Eq.  (3.61). 

J:K)-<7e=0  (3.61) 


Introducing  Eqs.  (3.59)  into  Eq.  (3.61)  yields 

rd2u  |  du    k2u  _  C23-C13c,0, 
dr2     dr        r  C 


(3.62) 
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+  cn-cu  a^T  +  c22-c,2 a  AT  +  c23-c13 a  AT 

C22  C22  C22 


where  the  anisotropy  factor  k  is  given  by 

*2  =7LL-  (3.63) 
c22 

The  general  solution  of  Eq.  (3.62)  may  be  written  in  the  form 

u  =  D/ +D2rk +Al£'°'r  +  A2r  (3.64) 

where  D\  and  D2  are  integration  constants,  and  A]  and  A2  are  constants  based  upon  the 
coil's  mechanical  and  thermal  properties. 

1     C  -C 

1    \-e  c22  (3-65) 

1  (c,2-c„)«eAr+(c22-c12)«rAr+(c23-c13)«7Ar 

For  an  isotropic  material  (k  =1),  the  right  hand  side  of  Eq.  (3.62)  vanishes,  and  hence,  A] 
and  A2  become  zero. 

The  solution  for  the  displacement  introduced  through  Eq.  (3.59)  yields  the  state  of 
stress  and  strain  as  given  in  Eqs.  (3.66)  -  (3.70),  together  with  Eq.  (3.60). 

ee  =  D/-'  +  D2r-*-'  +  \e?  +A2-  aeAT  (3 M) 

er  =  Dy~x  -  D2kr-k  X  +  A,ef  +A2-  arAT  (3.67) 

a,  ={Cu+kCi2)D/-i+{Cu-kCi2)D2r-k-]  (3.68) 
+(C„  +  C12)A,ef  +  (C„  +  C^A,  +  CI3< 


,/0/ 

'-z 


-CnaeAT  -  CnarAT  -  C„azAT 


ar  =  (C12  +  kC22)D/~x  +  (C12  - kC22)D2r-k~l  (3.69) 

+(C]2  +  C22)A£f  +(C12  +  C22)A,  +  C23< 
-C12aeAr  -  C22arAT  -  C23azAT 
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Oz  =  {Cu  +  kC2i)D/-<+(Q3-kC23)D2r^  (3.70) 

+(c13  +  c2,)a^T + (cI3  +  c23)^  +  c33£;<» 

-C13aeAr  -  C23arAT  -  C33azAT 

These  equations  together  with  the  value  of  the  integration  constants  D\  and  D2  and  the 
value  of  the  constant  axial  strain  ef,  determine  the  stress  and  strain  in  the  coil.  The 
equations  for  the  integration  constants  and  the  axial  strain  are  given  in  Appendix  A. 


CHAPTER  4 
GREEN'S  FUNCTIONS  SOLUTION 

Introduction 

A  full  axisymmetric  stress  analysis  in  terms  of  Green's  functions  for  an  elastic 
isotropic  coil  of  a  superconducting  magnet  is  presented.  This  analysis  is  applicable 
throughout  the  coil  and  includes  the  shear  stress  in  addition  to  the  normal  stresses.  As  the 
size  and  field  strength  of  high  field  magnets  increase,  the  knowledge  of  stress  distribution, 
especially  shear  stress,  beyond  the  midplane  of  the  coil  becomes  increasingly  important.  A 
limited  number  of  publications  on  the  approximation  of  the  three-dimensional  problem  are 
available  in  the  literature.22"25  Some  solutions  are  obtained  by  neglecting  shear  throughout 
the  coil.22' 23  Other  solutions  are  provided  based  upon  numerical  techniques24  or  power 
series  expansion  of  the  fields  and  displacements.25  The  direct  analytical  solutions  are 
formulated  for  infinite  or  semi-infinite  domains  that  are  not  appropriate  for  a  finite  domain 
such  as  a  magnet.26'27  It  is  desirable  to  develop  a  full  stress  analysis  for  a  superconducting 
magnet  irrespective  of  the  type  of  field  or  its  distribution  within  the  coil.  As  a  first 
approach  a  coil  with  isotropic  material  properties  is  assumed. 

In  the  following  sections  stresses  are  formulated  using  stress  functions.  Stress 
functions  are  solved  in  terms  of  Green's  functions  by  exploiting  finite  integral  transforms, 
e.g.,  finite  Hankel  transform. 

Fundamental  Equations  for  the  Stress  Functions 

Consider  an  elastic  isotropic  coil  with  inside  radius  of  a,  outside  radius  of  b,  and 
length  of  2L  as  shown  in  Fig.  4.1.  For  an  axisymmetric  distribution  of  Lorentz  force 
density,  X(r,z),  the  equilibrium  equation  may  be  written  as 

25 


26 


v.a+x  =  o 


(4.1) 


where  (7  is  the  stress  tensor. 

The  constitutive  equation  relating  the  stress  tensor  to  the  strain  tensor  8  is  given  by 

a  =  X  fr(£)I  +  2\iZ  (4.2) 

where  A  and  fx  are  (Lame's)  elastic  coefficients  and  tr(£)  is  the  trace  of  the  strain  tensor. 


Figure  4. 1  Schematic  diagram  representing  one  coil  of  a  magnet. 

The  displacement  vector  U(r.t)  is  related  to  the  strain  tensor  as  described  by  Eq.  (4.3). 

£  =  i[Vu  +  VTu]  (4.3) 
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Incorporating  Eq.  (4.3)  into  Eq.  (4.2)  yields  an  equation  for  the  stress  tensor  in  terms  of 
displacement  vector. 

a  =  X  (V.u)I  +  fi[Vu  +  VTu]  (4.4) 

Introducing  Eq.  (4.4)  into  Eq.  (4.1)  results  in  a  partial  differential  equation  for 
displacement  vector. 

(A+/x)V(V.u)  +  ^V2u  +  X  =  0  (4.5) 

From  the  Helmholtz  theorem,  any  vector  satisfying  Eq.  (4.5)  may  be  resolved  into 
a  sum  of  a  gradient  and  a  curl 

u  =  V0  +  VxA  (4.6) 

where  0<r.?)  is  a  scalar  potential  and  Ao-,*)  is  a  vector  potential  such  that  V.A  =  0. 
Incorporating  the  displacement  vector  from  Eq.  (4.6)  into  Eq.  (4.5)  yields  an  equation  in 
terms  of  potential  functions  <p  and  A. 

(A  +  2jU)V(V>)  +  ^Vx(V2A)  +  X  =  0  (4.7) 

It  is  admissible26  to  express  the  independent  potential  functions  0  and  A  as 

A  =  aVxvF  0  =  0V.*F  (4.8) 

where  a  and  p  are  arbitrary  constants  and  components  of  the  vector  *F  are  the  stress 
functions.  Introducing  Eq.  (4.8)  into  Eq.  (4.7)  leads  to  a  partial  differential  equation  for 
vector  *F. 

P(X  +  2//)V4lF  +  [pa  +2[i)  +  na]V  x  [V  x  (V2lF)]  +  X  =  0  (4.9) 

In  order  to  simplify  Eq.  (4.9),  we  may  choose  arbitrary  constants  a  and  P  as  — -  and 

- — — ,  respectively.  Thus,  Eq.  (4.9)  reduces  to 
A  +  2./I 
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v4XF  +  x  =  o 


(4.10) 


In  cylindrical  coordinates  ( r,  6,  z),  Laplacian  of  vector  *F  is  given  by 


V2XF  = 


v/2  uj     Y,     2  d*Fe 
'     r2     r2  dd 


e.  + 


V2^ 


r2     r2  96 


ee  +  V2lF2ez  (4.11) 


where  the  Laplacian  ( V2)  of  a  scalar  is  represented  by 


v2=il  ii.  1  d%  <L 

dr2  +  rdr  +  r2  d62  +  dz2 


Applying  Eq.  (4.11)  to  Eq.  (4.10)  results  in  three  partial  differential  equations  for  stress 
functions.* 


i  V 


V  "X 
r  J 


r  J 


Yc  - 


4  d2x¥r  4 


r4  dd2  r2 


4  d2Va  4 


e  r4 


<?02  r2 


2 


50 


+  X  =0 


(4.12) 


+  xe  =  o 


V4«P:  +  XZ  =0 

The  geometry  is  axisymmetric  in  nature  and  for  an  axisymmetric  system,  the  partial 

derivative  with  respect  to  tangential  direction  is  zero,  —  =  0.   Thus,  the  three  partial 

d6 

differential  equations  in  Eq.  (4.12)  will  reduce  to  three  uncoupled  partial  differential 
equations  given  by  Eqs.  (4.13)  -  (4.15). 


1  \2 

2  1 


V2- 


r2j 


•P.  +  X  m  0 


r  r 


(4.13) 


(  ,   1 V 

'v2  — U  f0+xe  =  o 


(4.14) 


See  Appendix  B  for  details  of  the  derivation. 
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V4f2+Xz=0  (4.15) 

The  driving  body  force  in  a  magnet  is  due  to  the  Lorentz  force  density.  The  Lorentz 
force  is  a  function  of  r  and  z,  and  is  related  to  the  magnetic  field  and  current  density  by 

X  =  JxB  (4.16) 

where  J  and  B  are  the  current  density  and  the  magnetic  field  vectors,  respectively.  For  an 
axisymmetric  magnet  J  =  Jgee  and  B  =  Brer  +  Bzez.  The  vector  product  of  the  J  and  B 

leads  to 

Xr  =  JeBz  Xe=0  Xz=-JeBr.  (4.17) 

Thus,  in  the  absence  of  a  tangential  magnetic  body  force,  Eq.  (4.14)  yields  *Fg  =  0,  which 
results  in  tangential  displacement  ug  to  be  zero.  Solutions  to  Eqs.  (4.13)  and  (4.15)  will 
determine  the  stress  functions  and  hence  the  stresses.  In  order  to  solve  these  equations,  we 
need  to  exploit  the  finite  Hankel  transform. 

Finite  Hankel  Transform 

In  contrast  to  the  Hankel  transform  that  is  available  in  the  literature  corresponding  to 
an  infinite  domain,  the  finite  Hankel  transform  is  a  rare  commodity.28- 29  The  ones  that  are 
available  are  mainly  used  for  an  interval  starting  at  the  origin.  Hence,  they  carry  only  the 
Bessel  function  of  the  first  kind.  For  the  case  of  a  coil  (a  hollow  cylinder),  the  radial 
interval  starts  at  a  point  greater  than  zero,  thus,  a  suitable  finite  Hankel  transform  must 
contain  the  Bessel  function  of  the  first  kind  and  the  Bessel  function  of  the  second  kind.  A 
detailed  derivation  of  the  finite  Hankel  transform  for  an  interval  appropriate  for  a  coil  is 
carried  out  in  this  section. 

Derivation 

A  function  //(a),  as  expressed  in  Eq.  (4.18),  defines  an  integral  transform  of  an 
arbitrary  function  /<*>  by  a  kernel  K«x,x). 
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Ifw  =  j  f(*)Kia,x)dx  (4.18) 
o 

The  method  of  integral  transform  may  under  some  circumstances  be  applied  to  the  solution 
of  boundary  value  and  initial  value  problems  in  mathematical  physics.  In  problems  in 
which  one  of  the  independent  variables,  say  x,  has  a  range  of  (0,°°)  the  use  of  an  integral 

transform  of  the  type  j  fu)K(a,x)dx  will  reduce  a  partial  differential  equation  in  n 

o 

independent  variables  to  n-l  independent  variables  thus,  reducing  the  complexity  of  the 
problem.  In  some  instances  successive  operations  of  this  type  will  ultimately  reduce  the 
partial  differential  equation  into  an  ordinary  differential  equation.  It  is  possible  to  extend 
the  method  of  integral  transforms  to  a  finite  interval  in  which  the  field  of  variation  of  the 
independent  variable  is  now  [a,b]  where  both  a  and  b  are  real  and  finite. 

A  Fourier  Bessel  kernel  Kn((,,x)  based  upon  a  solution  of  a  Bessel  equation, 
1      (  n2\ 

y  =  0,  with  homogeneous  Dirichlet  conditions  in  the  finite  interval 
[a,  b]  where  a  is  greater  than  zero,  is  expressed  by 

*s,<f,.*)  =  [4,  /.(^.(O)].  (4.19) 

Here,  £(.  is  a  root  of  the  transcendental  equation  given  by  Eq.  (4.20). 

Jn((,a)Yn^b) -  J^cpY^a)  =  0  (4.20) 

The  kernel  Kn((.,x)  is  orthogonal*  in  the  interval  [a,b]  with  respect  to  a  weighting  function 
rw  =  x .  Thus,  we  may  write 

J^.^«,.^-(*     f  (4-21) 


y  +—y  + 

r 


C-n- 

r  J 


See  Eq.  (C.18)  of  Appendix  C. 
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where  both  £.  and  are  roots  of  Eq.  (4.20)  and  ||£„(f,,*)|  is  the  norm  of  the  Kn(t,,x)  that 
is  expressed*  by  Eq.  (4.22). 


V2 


l- 


(4.22) 


From  the  theory  of  Bessel  functions,  if  satisfies  the  Dirichlet  conditions  in  the  interval 
[a,b],  then  the  series 


converges  to  the  sum  of  —  [/(*  +  o)  +  /<*  -  o>].  Thus,  if  fw  is  continuous  at  the  point  x  the 

2 


series  converges  to  /<*). 


(4.23) 


In  Eq.  (4.23)  A«;.)  is  an  arbitrary  function  of  A^^.x)  is  the  Fourier  Bessel  kernel,  and 
the  summation  is  taken  over  all  the  positive  roots  of  £, . 

In  order  to  obtain  an  expression  for  the  A(f.)  in  Eq.  (4.23)  we  make  use  of  the 
orthogonality  property  of  Kn(crx).  Multiplying  both  sides  of  the  Eq.  (4.23)  by  the  term 
xKni(Jtx),  integrating  over  the  interval  [a,b],  and  assuming  that  term  by  term  integration  is 

permissible,  yields 


D  D 

]xf(x)Kn(srx)dx  =  ^A(i;i)jxKn(C.,x)Kn((J,x)dx. 


(4.24) 


All  terms  on  the  right  hand  side  of  Eq.  (4.24)  vanish  except  the  term  corresponding  to 
i  =  j.  Thus,  Eq.  (4.24)  is  reduced  to  Eq.  (4.25). 


See  Eq.  (C.38)  of  Appendix  C. 


32 


jxf(x)Kn(zrx)dx  =  A^K^.x-f  (4.25) 

a 

Eq.  (4.25)  is  the  foundation  for  the  finite  Hankel  transfer  theorem.  If  /(*)  satisfies 

Dirichlet  conditions  in  an  interval  [a,b],  then  its  finite  Hankel  transform  of  order  n,  denoted 
by  /({,),  is  defined  by 


b 

%[f^]  =  /(£,)  =  jxfu)Kni(rx)dx  (4.26) 
a 

where     is  a  root  of  the  Eq.  (4.20),  9t„  [/(*>]  is  a  linear  functional  operator,  and  J^.x) 

is  the  Fourier  Bessel  kernel  given  by  Eq.  (4.19).  By  using  Eqs.(4.25)  and  (4.26)  we  may 
write  the  arbitrary  function,  A<c(),  in  terms  of  the  finite  Hankel  transform  of  /<*). 

A(;i)  =  \\K^,x)l~2f(Q  (4.27) 

Incorporating  Eqs.  (4.27)  and  (4.22)  into  Eq.  (4.23)  yields  the  following  inversion 
theorem.  If  /(*>  satisfies  the  Dirichlet  conditions  in  an  interval  [a,b]  and  if  its  finite  Hankel 
transform  is  fa.  >  then  at  each  point  in  that  interval  fix)  may  be  written  as 


/«  =  9C[fa>]  -  TX-^-MM  (4-28) 

where  the  summation  is  extended  over  all  the  positive  roots  of  . 

Finite  Hankel  Transforms  of  the  Derivatives 

Suppose  that  fa.)  is  the  finite  Hankel  transform  of  /<*>.   The  finite  Hankel 

transform  of  the  one-dimensional  Laplacian  of  /u>  in  cylindrical  coordinates, 
f     1  df 

V2/  =  —4  +  — i-,  is  then  expressed  by  Eq.  (4.29). 
dx      x  dx 


5t  [£l+L4L]=]x 

dx2     x  dx  J 


<*7     1  4f  , 
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Integration  by  parts  of  the  right  hand  side  of  the  Eq.  (4.29)  gives 

b 


d'f  ,  1  df 


dx      x  dx 


df 


+  ——  Kni(..x)dx  =  x-^-Kn((.,x) 


dx 


-xf 


dK. 


dx 


iJ  dx 


X) 


dx 


dx. 


(4.30) 


The  kernel  Ka%,x)  in  Eq.  (4.30)  is  derived  from  Bessel  functions  of  the  first  and 
second  kind,  therefore,  it  satisfies  the  Bessel  equation. 


dx 


dK, 


dx 


]  +  (£V-n2)tfn(c„,>  =  0 


(4.31) 


Introducing  Eq.  (4.31)  into  Eq.  (4.30)  yields 


d2f  ,  1  df 


dx      x  dx 


df 


+  — —  Kn(;.,x)dx  =  x— Kn((.,x) 


dx 

h 


-xf 


x) 


dx 


b     f  2  \ 


Kn(trx)dx 


(4.32) 


or 


d2f    Idf  n7 


+  ;  r/  K„((l.x)dx  =  X  —  Kni(rx) 

dx      x  dx  x 


dx 


-xf 


dK 


dx 

b 

-  Q  J xf(x)KH(crx)dx. 


(4.33) 


The  first  term  on  the  right  hand  side  of  Eq.  (4.33)  vanishes  at  both  upper  and  lower  limits. 
The  integral  part  of  the  last  term  is  the  finite  Hankel  transform  of  /(*>.  Hence,  Eq.  (4.33) 
takes  the  following  form. 
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Now,  it  is  shown  that" 


Xf(x) 


dK 


dx 


(4.35) 


The  substitution  of  Eq.  (4.35)  into  Eq.  (4.34)  yields  an  expression  for  the  finite  Hankel 
transforms  of  the  derivatives  of  /<*>. 


ax      x  ax  x 


K 


fW-~f-J-fla) 


(4.36) 


Green's  Functions 

The  partial  differential  equation  for  radial  and  axial  stress  functions  are  given  by 
Eqs.  (4.13)  and  (4.15).  Solutions  to  these  equations  may  be  obtained  by  variety  of 
methods.  The  present  work  provides  the  solutions  in  the  form  of  Green's  functions.  The 
usage  of  Green's  functions  allows  us  to  obtain  a  solution  irrespective  of  the  field's  type  or 
its  distribution  within  the  coil.  Moreover,  the  same  Green's  functions  may  be  used  in  order 
to  obtain  a  solution  for  the  thermal  stresses. 

In  the  next  two  sections,  solutions  to  Eqs.  (4.13)  and  (4.15),  with  Dirichlet 
homogeneous  conditions  are  obtained  in  the  form  of  radial  and  axial  Green's  functions, 
respectively. 

Radial  Green's  Function 

In  order  to  find  a  solution  for  Eq.  (4.13),  two  finite  integral  transforms,  the  finite 
Hankel  transform  in  r  and  the  finite  Fourier  cosine  transform  in  z,  are  exploited.  Applying 
the  finite  Hankel  transform  (of  order  one)  to  each  term  of  the  Eq.  (4. 13)  yields 

^[[v'-ijyj-ft.fx,].  (4.37) 
*  See  Eqs.  (C.35)  and  (C.36)  of  Appendix  C. 
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By  using  the  property  of  the  finite  Hankel  transform  (for  derivatives),  Eq.  (4.37)  reduces 
to 


22  \2 
fc2 


(4.38) 


where  *Fr((u.z)  and  Tr\(u,t)  are  the  finite  Hankel  transforms  of  the  radial  stress  function 
and  radial  body  force,  respectively. 


(4.39) 


Xr  c{s . ,  z)  =  J  rXr(r,  Z)Kt  (fu .  r)dr 


(4.40) 


The  Fourier  Bessel  kernel  jr,(fK,r)  in  Eqs.  (4.39)  and  (4.40)  is  expressed  by 


where  £u  is  a  root  of  the  corresponding  transcendental  equation. 


(4.41) 


(4.42) 


An  additional  transform  (in  the  axial  direction)  is  imperative  for  solving  Eq.  (4.38). 
Considering  the  radial  body  force  is  an  even  function  of  z,  an  appropriate  transform  is  the 
finite  Fourier  cosine  transform.  By  introducing  the  finite  Fourier  cosine  transform 


dz2 


(4.43) 


the  ordinary  differential  equation,  Eq.  (4.38),  is  converted  into  an  algebraic  equation 
expressed  by 


-2       n  K 

"feu  7T 


Yr(tu.H)  =  -Xri(u,n) 


(4.44) 
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where  n  is  an  integer  number.  Functions  fr and  Xr(C„.«)  are  the  finite  Fourier  cosine 
transforms  of  Y^.t)  and  Xrau,t),  respectively. 

  z. 

!Fr«^>«  jWr{(u.z)oos^dz  (4.45) 

_  L 

Xr(c1(,«)=  JXr(fh,Z)Cos^-^Jz  (4.46) 
From  Eq.  (4.44),  ?Fr({^o  may  be  written  in  terms  of  X^.n). 

^u.-)-r      "L'     ,1^,,)  (4.47) 

The  inverse  finite  Fourier  cosine  transform  of  ?Fr ({;„,»)  is  defined  by  Eq.  (4.48). 


^.0»3^T^,(fw.")]  =  ^^r(C».o)  +  7]£^--)<»^  (4.48) 

^  „=i 


Incorporating  Eq.  (4.47)  into  Eq.  (4.48)  yields 


*r«r*  =  ~^Xr^-L'±  1  2         Xr,CB„)COS^.  (4.49) 


The  inverse  finite  Hankel  transform  of  *Pr(iu.t)  is  given  by  Eq.  (4.50),  where  Kj({H.r)  is 
the  same  Fourier  Bessel  kernel  expressed  by  Eq.  (4.41). 

2  t2 


= 9C[ *X-J = f-I  ,2  M  Y    yr«;.^(i..»  (4.50) 
2  ,.=1  y,  (fua)-7,  (fuw 


The  substitution  of  fr(?„..Z)  from  Eq.  (4.49)  into  Eq.  (4.50)  yields  an  equation  for  the 
radial  stress  function. 
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CiVi2(c„") 


2  w  y,  (f„.fc) 


1  V 


(4.51) 


n(L2^+n2^2) 


j-Xf.(fl(.ii)Cos 


Introducing  Eqs.  (4.40)  and  (4.46)  into  Eq.  (4.51)  gives  the  final  solution  of  the  radial 
stress  function  *Fr(r,z)  in  terms  of  the  radial  body  force  Xr(r,Z). 


1    /=1  [Z^bl,  -ta 


(4.52) 


n=\ 


-cos- 


hkz 


L  b  , 

J  J  r'Xr(r',      (fu ,  r)  cos  dr'dz' 


-La 


2  2  A.,^..!-), 


By  interchanging  integrals  with  summations,  we  can  write  Eq.  (4.52)  in  the  form 


2LR 


(4.53) 


nitz  nnz 

-cos  cos  

L  L 


>XAr\z')dr'dz' 


or 


*Fr(r,z)=  Xr(r',z')G  r(.r,r',z,z')dr'dz' 

J-LJa 


(4.54) 


where  Gr(r,r\z,z')  is  the  radial  Green's  function  defined  by  Eq.  (4.55). 


Gr(r,r',z,z')  =  -^ 


-  -  j  ~2  r  A,  (f , . ,  r) ^  (f , . ,  r') 


,=i  I  2  y,  (f„.a)-y,  <?„.« 


i 


+ 


cos  cos 


2LQ       tt(L2Q  +  n2x2f        L  L 


(4.55) 


Axial  Green's  Function 
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The  axial  Green's  function  is  obtained  by  solving  the  partial  differential  equation  for 
the  axial  stress  function,  Eq.  (4.15).  In  order  to  solve  Eq.  (4.15),  the  same  finite  integral 
transform  (finite  Hankel  transform)  with  a  different  order  in  r  is  used.  However,  since  the 
axial  body  force  is  an  even  function  of  z,  the  finite  Fourier  sine  transform  (instead  of  finite 
Fourier  cosine  transform)  is  utilized  in  the  axial  direction. 

Applying  the  finite  Hankel  transform  of  order  zero  to  each  term  of  the  Eq.  (4.15) 

yields 


By  using  the  property  of  the  finite  Hankel  transform  (for  derivatives),  Eq.  (4.56)  reduces 
to 


where  !Fz<cw.z)  and  Xz(f0j,z)  are  the  finite  Hankel  transform  of  the  axial  stress  function  and 
axial  body  force  given  by  Eqs.  (4.58)  and  (4.59),  respectively. 


(4.56) 


( 


d2  V  — 


(4.57) 


b 


(4.58) 


a 


b 


(4.59) 


a 


Here,  K0(torr)  and  £0(  correspond  to  zero  order  transformation. 


(4.60) 


(4.61) 
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By  taking  the  finite  Fourier  sine  transform 


(  #\ 


3S[  -a+tr  ?,«;.«)]— 3fR«;.o] 


dz2 


(4.62) 


the  ordinary  differential  equation,  Eq.  (4.57),  is  converted  into  an  algebraic  equation 
expressed  by 


(  „2tt2\ 


(4.63) 


where  n  is  an  integer  number,  and  *Ft{Cm,*)  and  X^,*)  are  the  finite  Fourier  sine 
transforms  of  the  ¥ttfv,i)  and  Xt(CK.z). 


(4.64) 


-L 


From  Eq.  (4.63)  we  can  write  *Pz(f0l.«)  in  terms  of  X. ({„,»). 


(4.65) 


-L4 


^•")=(^+«v) 


2  X,(f«.») 


(4.66) 


The  inverse  finite  Fourier  sine  transform  and  finite  Hankel  transform  are  defined  by 
Eqs.  (4.67)  and  (4.68),  respectively. 


(4.67) 


^=9g[^»z>]^£  /2  ay°2(C;;>  (4.68) 

2  J0(C0M)-J0(S0ih) 
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Incorporating  Eqs.  (4.66)  and  (4.67)  into  Eq.  (4.68)  yields  an  expression  for  the  axial 
stress  function. 


--XI 

1=1  n=l 


2^ 


(L2a+»v) 

K1  I1 


(4.69) 


nnz 


Jo  <fo,«>  -  -A)  tV) 


Xr(f0i,n)sin— -K0(corr) 


Introducing  Eqs.  (4.59)  and  (4.65)  into  Eq.  (4.69)  gives  the  final  solution  to  the  axial 
stress  function  W^ir.i)  in  terms  of  the  axial  body  force  Xr(r,Z). 


-Ill-, 


(4.70) 


Lbr 

w 


-La 


nnz 


r'K0(Z0i,r')sm——Xz  (r-,z-)dr'dz 


r)sm  - 


nnz 


By  interchanging  integrals  with  summations,  Eq.  (4.70)  may  be  written  in  the  form 


L  b 


-JJXI 

-La     '='  «=1 


7T 


(S0iJ0  (C0,a> 


2  (L2C+«Vf  *«^-*«; 


(4.71) 


„  n7Ez   .  nnz 

r  K^.r'yK^.rKm—j-sm  — 


>X<r:Z')dr'dz' 


or 


(4.72) 


where  Gz(r,r-.z.z')  is  the  axial  Green's  function  defined  by  Eq.  (4.73). 


G!(ry.«')  =  --JJ 


2  ,=i  «=i 


(L2£2  +  n  V)      <fo,a> "  -A) 

„        .  nnz'  .  nnz 
r  K0(c0, . n K0 (f0i . rjsin  — —  sin  — 


(4.73) 
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Boundary  Conditions 


For  traction  free  boundary  conditions,  the  radial  and  shear  stresses  are  zero  at  the 
inside  radius,  r  =  a,  and  outside  radius,  r  =  b,  and  the  axial  and  shear  stresses  are  zero  at 
the  ends  of  the  coil,  z  =  ±L.  Since  the  traction  free  boundary  conditions  were  not 
considered  for  the  solution  of  the  stress  functions,  the  expression  for  the  stresses  obtained 
from  Wr(r,Z)  and  *Ft(r,z)  will  not  satisfy  these  boundary  conditions.  Thus,  radial  and  shear 
stresses  impose  new  non  zero  forcing  functions  of  z  at  radial  boundaries,  and  similarly 
axial  and  shear  stresses  assert  new  forcing  functions  of  r  (in  general  non  zero)  at  axial 
boundaries.  To  alleviate  this  apparent  difficulty,  other  solutions  (complementary  solutions) 
for  stress  functions  must  be  obtained  in  order  to  reverse  the  effect  of  these  new  imposed 
conditions.  Since  the  partial  differential  equations  for  stress  functions  (Eqs.  (4.13)  and 
(4.15))  are  linear,  the  superposition  principle  is  applicable.  Hence,  the  combination  of  the 
complementary  solutions  and  the  Green's  functions  furnish  the  final  solutions  for  the  stress 
functions.  In  this  section  the  imposed  forcing  functions  are  derived,  and  in  the  succeeding 
section  the  complementary  solutions  are  developed. 

The  displacement  vector  is  related  to  the  vector  *F  from  Eqs.  (4.6)  and  (4.8). 


The  curl  of  the  curl  of  any  vector  may  be  expressed  by  the  gradient  of  the  divergence  of  that 
vector  minus  the  Laplacian. 


Incorporating  Eq.  (4.75)  into  Eq.  (4.74)  and  substituting  *F  by  +  fs(r.i)ea  into 

Eq.  (4.74),  yields  an  equation  for  the  displacement  vector  in  terms  of  the  stress  functions, 
«Pr(r.z)  and  fz<r,z). 


u  =  aVx(VxKF)  +  )8V(V.vF) 


(4.74) 


(4.75) 


u  =  -aV2(tPr(r,z)er+  «F,(r.Z)eI)  +  (a  +  )8)V[V.CPr(r,J,er  +  "F^eJ  (4.76) 
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Rewriting  Eq.  (4.76)  in  component  form  results  in  expressions  for  radial  and  axial 
displacements,  where  <jP  is  the  divergence  of  vector      given  by  Eq.  (4.78). 


u  =  -a\  V2  - 


r2J 


uz  = 


dr 

d(p 


-aV2f,+(a  +  /?)^- 
dz 


(4.77) 


1  d  ,  m  .  *IP, 


(4.78) 


The  stress  tensor  in  terms  of  the  displacement  vector  is  defined  by 

G  =  A(V.u)I  +  /x[Vu  +  VTu]. 


In  the  component  form  it  becomes 

<7  =  A 


CTe=A 


G.  =  A 


1  d  .  .  du7 

.  r  or  cte 

1     >  \  duz 

7</r  dz 

Id.  .  du7 

r  dr  dz 


dr 


_  u 
+  2yL  — 
r 


~  du. 
+  2fi — z- 


dz 


dur  du 


dz  dr 


(4.79) 


The  substitution  of  displacements  from  Eq.  (4.77)  into  radial  stress  from  Eq.  (4.79),  yields 


j     1  d 
<yr=-Xa-— 

r  dr 


.  \        r  J 

1  d  (  d<p\ 


r  dr 


r 

v  dr  ) 


-A«|(V^z)  +  A(a  +  ^^ 
dz  dz 


\dz  j 


dr 
d 


(4.80) 


1  ^ 


IP. 


+  2n(a  +  /3) 


d_(d£ 
dr\dr  > 
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The  terms  -— 
r  dr 


{r-7) 

r 

and 

<        r  ) 

and 


(  d  \ 

^dz  % 


,  respectively.  Hence,  substitution  of  these  terms  into  Eq.  (4.80)  gives 


cr  =  -A  aV' 


-2tia-r\  V2-- 

dr\        r  ) 


1 


fr  +  2/*(a+/3) 


r  dr 

dr 


1  d  (Jq>\  |  <?> 


\  dr  j 


dz2 


(4.81) 


\dr  j 


\d_ 
r  dr 


dVz  . 
dz 


1  d  (  d(p\  d2q> 


but  — —(r*Fr)  +  —^L  is  <j9  and  — —  r^  is  the  Laplacian  of  (p.   Thus,  Eq. 

r  dry  dr  )  dz 


(4.81)  becomes 


ar  =  A  pV2(p  -  2/ia 


dry        r"  J 


fr  +  2/i(a  +  /J) 


d2q> 
dr2  ' 


(4.82) 


The  substitution  of  a  =  — -  and  B  =  — - —  in  Eq.  (4.82),  yields 

//  A  +  2/J. 


or  = 


r    A  +  2^ 


V>  +  2 


d_ 
dr 


f*4) 

V        r  ) 


^  +  2(-l  +  -iL_)^f 
A+2// <?r2 


(4.83) 


Lame's  elastic  coefficients,  A  and  /i ,  are  related  to  each  other  by 

2//v 


A  =- 


l-2v 


(4.84) 


where  vis  the  Poisson's  ratio.  Introducing  Eq.  (4.84)  into  Eq.  (4.83)  gives  the  radial 
stress  in  terms  of  stress  functions. 


a  =2 


I  r2) 


dr 


1 


1- v 


vV2- 


i_2> 

dr  j 


(4.85) 
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In  a  similar  fashion,*  other  components  of  stress  are  determined  in  terms  of  stress 
functions. 


°e=- 


(  1  "\ 
V        r  J 


1 


r    1- V 


vV2- 


V 


ld_ 

r  dr  > 


(4.86) 


vV2- 


<?z2 


r2J 


1  d2cp 
1  -  V  drdz 


V- 


Stress  functions  fr(r,Z)  and  f  (r,t)  are  defined  by  Xr(r',z-)Gr(r,r',z,z')dr'dz/ 


and  J  J  Xz(r-,z')Gz(r,r',z,z')dr'dz'  where  Gr(r,r-,z,z')  and  Gt(r,r',«,t')  are  radial  and  axial 


Green's  functions  given  by  Eqs.  (4.55)  and  (4.73),  respectively.  Based  upon  Eqs.  (4.55) 
and  (4.73)  and  by  introducing  new  functions  rr(fB,«)  and  r  <{„,«),  the  stress  functions 


may  be  written  as 


Yr  <r'  z>  =  X  X  r^ (f«  •  n)^i  Cf«  • r) C0S 


i=l  n=0 


(4.87) 


1=1  n=l  ^ 


(4.88) 


where  Tr{(v.n)  and  rz<cw,«)  are  given  by  Eqs.  (4.89)  and  (4.90). 


rt(f0(,»)=- 


T 


Coi-A)  <£p,a> 


(L2Co,-  +  n  V)  4>  «o,«> -  Jo<V> 
f   f  r'Xz(r-,z')  ^0(f0i,r-)sin^-^r'rfz' 


(4.89) 


*  See  Appendix  D  for  details. 
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K1 


rr(f„.o)=- 


2  JU^-JU^tl+rtK1)1 

\L  ^rXir^K^.ncos—dr'dz' \ 

J-LJa  L  J 

—TT\   \  r'Xr(r:z-)K^hy)dr'dz' 

.L.  *-LJa 


2  7,  (?„*.)  2L£1(. 


(4.90) 


The  substitution  of  Eqs.  (4.88)  and  (4.87)  into  Eq.  (4.85)  gives 


*=*iXr,«,.>|(v>-^ 


i=l  n=0 
1 


co  OO 


vV2- 


dr 


vV2- 


j 


r  dr\_ 


nnz 

rKj(fK,r)COS— — 


tf0(C0i.r)Sin— - 


(4.91) 


The  Fourier  Bessel  kernel  Jr,(CB.r)  is  obtained  from  Bessel  functions  of  order  one. 
Hence,  it  satisfies  the  Bessel's  equation  which  translates  to 


f 

1  " 

r  dr 

r 

< 

dr) 

r2_ 

K]((li,r)  =  -£iKl(iirr). 


(4.92) 


Also  from  basic  identities  for  Bessel  functions,  we  obtain 


feu'' 


or  in  terms  of  JT,(Cu,r) 


1 


(4.93) 


(4.94) 
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Likewise  for  K0^g.,r)  we  may  write 


r  dr 


f  — 


(4.95) 


Using  results  of  Eqs.  (4.92),  (4.94)  and  (4.95)  and  by  introducing  new  kernels  K0}{(  r) 


and  Klo((orr),  we  can  reduce  Eq.  (4.91)  to 


oo  oo 


n) 


1=1  n=0 


-a 


2_2  \ 


n  K 


I 


\ 


r  J 


cos- 


nnz 


(4.96) 


1  -  V  i=\  n=0 


V 


zV2 


~2     W  7T  «- 


nnz 
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V  r  J 


cos- 


ttTTZ 


+  t— ;IXr«<f.-") 

1       V  i=l  n=0 


v 


V 


-C2  - 

bo; 


2_2  \ 


2    n  7T 

Oi 


£0(f0i.r)COS 


L) 


1 


Co,  Co1^0(fo,-r»~_^lo(?o,-'-> 

V  r  J 


cos- 


where  15roi(fu.r)  and  A"10(c0i,r)  are  given  by  Eq.  (4.97)  and  (4.98). 


1  1  d 

Km^r*  =  [JniWYx«ub)- Jx(tub)Y0«nr)]  =  —Kl«u.r)  +  ——Ki{(v,r)  (4.97) 


1  <? 


(4.98) 


Eq.  (4.96)  may  be  simplified  further  by  combining  the  summations  and  canceling  the 
common  terms. 
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z=l  n=0 


T 


v^r-fi(i-v) 


(4.99) 


«7T 


-rr(f|»--)Cli 

+rr  (?„,«) 


2  2 

n  7T 


(2-v)+fi(l-v) 


2^1(l-v)+fi(l-2v) 


K^irr)\  COS 


Likewise*  we  may  write  the  other  components  of  stress  directly  in  terms  of  the 
Fourier  Bessel  kernels  and  trigonometric  functions. 


OO  OO 


n  K 


rK0aorr) 


(4.100) 


H7T 


-rr(cu.») 


-72-+ C- 


2^(l-v)+C,-(l-2v) 


^(^.D^COS 


n;rz 


CT-||T^i-r-,c-'f 


2  2 

n  7T 


(l-v)  +  Co2,(2-v) 


K^0i,r)  (4.101) 


2  2 

n  K 


(l-v)-<. 


^(^.D^COS 


n;rz 


i=\  n  =  \ 


V^-Co2,(l-v) 


T  r<<:"'n) 


2  2 


-(l-v)-< 


„         .  nnz 

K^,r)\SlX\— — 


(4.102) 


The  above  expressions  for  stresses  satisfy  the  field  equation  but  not  the  boundary 


*  See  Appendix  E  for  details. 
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conditions.  At  the  inner  radius,  by  the  use  of  K]((l.,a)  =  K0aol.a)  =  0,  and  from  Eqs. 
(4.103)  and  (4.104),  the  radial  and  shear  stresses  are  determined. 


10<fopa> 


(4.103) 


-rr(fB.»)Cii 


2  2 

n7r 


(2-v)'+fi(l-v) 


i=l  „=i  I    i  V 


v^-fi(l-v) 


„  nnz 
Ari0<fw.a)Sin— 


(4.104) 


From  Eqs.  (4.93)  and  (4.97)  we  may  write 


and 


Aio<?o,-a>-~7  


Co,  A- 


However,  AT,(fu,a)  is  zero  and  —  if^.r) 


<9r 


and  —  K0«„,r) 


may  be  expressed*  by 


 li£^  ^  pi£o^>  ^  reSpectively.  Thus,  Eqs.  (4.103)  and  (4.104)  are  reduced 


to 


w5o«a  1- v  I  L  70(f0,fl) 


(4.105) 


rr(fw,«) 


2  2 


(2-v)+e(i-v) 


cos- 


MtZ 


*  See  Eq.  (C.37)  of  Appendix  C. 
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1=1  «=i 


2  2 

«  7T        >-2  /i  \ 
 fi(l"  V) 


70(f0ife)  2     1     .  «7Tz 

-y— 21  sin  

70(f0ia)  sal-v  L 


By  introducing  new  functions  p,oo  and  £>2(«)  as 


(4.106) 


i=i 


_2  1_ 

7ra2  1-  v 


-r,(fw,n) 


2  2 


(2-v)+g(l-v) 


(4.107) 


n) 


•A><fo,<" 


(4.108) 


we  may  simplify  Eqs.  (4.105)  and  (4.106)  to  Eqs.  (4.109)  and  (4.1 10),  respectively. 


Or(a.z)  =  ^£>,(«)COS 


n=0 


(4.109) 


Xnnz 
p2(«)Sin — 


(4.110) 


At  outer  radius,  r  =  b ,  with  use  of  the  expressions* 


<9r 


,-=/> 


_2_ 
7T& 


and 


r=b 


_2_ 
Kb 


Eqs.  (4.105)  and  (4.106)  may  be  written  as 


See  Eq.  (C.37)  of  Appendix  C. 
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2  1 


+*rr(f„,«) 


«7T 


2  2 

n  7T 


(4.111) 


(2-v)+fl(l-v) 


•cos- 


MlZ 


oo  oo 


',*) = y.y. 
i=i  n=i 


1  2 


2  2 


V^~fi(l-V) 


•sin- 


(4.112) 


or 


#?3(«)COS  — 

n=0  ^ 


(4.113) 


v1 ,       .  nnz 


(4.114) 


where  p3w  and  p4(»)  are  given  by  Eqs.  (4.115)  and  (4.116). 


vi  2     1      _       nn    ,  _ 


2  2 


(2-v)+g(l-v) 


(4.115) 


1  2 
tl\-vnb\ 


2  2 

"  ^         >-2  /,  \ 
 C0,(l-  V) 


(4.116) 


ftJTZ 

Since  sin  vanishes  at  both  L  and  -L,  shear  stress  at  the  axial  boundaries  does 

L 

not  exert  a  forcing  function.  However,  axial  stress  is  a  function  of  cosine  and  does  not 
vanish  at  boundaries.  The  substitution  of  z  =  ±L  into  axial  stress,  Eq.  (4.101),  yields 


,=i  n=o 1    v  y  *- 


2  2 

n  n 


(l-v)  +  £(2-v) 


K0«0i.r)  (4.117) 


2  2 


(l-v)-vfi 


*01  (fw-O  I 
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or 


at(r,±L)  =  J [pj  (f0,  ^(f. ,  r)  +  ^6  <?,,         (fu ,  r)] 


where  p5(?0l)  and  p6(f„)  are  given  by  (4.119)  and  (4.120). 


(4.118) 


Ml 


2  2 

n  n 


(l-v)  +  Co2,(2-v) 


(4.119) 


n=0  1  "I 


2  2 

n  k 


(l-v)-v^ 


(4.120) 


Eqs.  (4.109),  (4.110),  (4.113)  and  (4.114)  together  with  Eq.  (4.118)  dictate  the 
forcing  functions  at  the  boundaries.  Complementary  solutions  for  stress  functions  must  be 
obtained  that  reverse  the  effect  of  these  forcing  functions.  Looking  at  the  equations  for 
stresses,  we  observe  that  all  of  the  components  of  stress  are  related  to  both  radial  and  axial 
stress  functions.  Hence,  either  one  of  the  stress  functions  may  be  used  to  reverse  the  effect 
of  the  forcing  functions.  Here  we  use  the  axial  stress  function.  Because  the  original  partial 
differential  equation  is  satisfied  by  the  axial  Green's  function  solution,  we  require  that  the 
complementary  solution  only  satisfies  the  homogeneous  part  of  the  partial  differential 
equation. 

Complementary  Solution  for  the  Axial  Stress  Function 

Consider  function  §<r.*)  as  a  complementary  solution  to  the  axial  stress  function. 
From  Eq.  (4.15),  £(r,z)  satisfies  the  homogeneous  partial  differential  equation  of  the  axial 
stress  function. 


VV;)  =  0 


(4.121) 
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From  Eqs.  (4.78),  (4.85),  and  (4.86),  radial,  axial,  and  shear  stresses  may  be  written  in 
terms  of  £(r.t). 


(7  = 


1  d 


"  l-vdz 


vV2- 


(4.122) 


1  d 


a  = 
;  l-vdz 


(2-  v)V 


<?2 


dz2 


(4.123) 


<7  = 


1  (9 


*  1-vdr 


(1- v)V2- 


^2 


£(r,z) 


(4.124) 


The  function  £(r,*)  must  reverse  the  effect  of  the  imposed  forcing  functions  by 
stresses  at  the  boundaries.  Thus,  from  Eqs.  (4.109),  (4.1 10),  (4.113),  (4.114),  (4.118), 
(4.122),  (4.123)  and  (4.124)  boundary  conditions  for  £ir,z)  may  be  expressed  by 


1  d 


l-vdz 


vV2- 


i_2> 

dr 


=  -Xp,<")COS^  (4.125) 

n=0  Li 


1  5 
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,<r,z) 


H  =  l 


(4.126) 


1  d 
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1  -  v  dz  \  dr 
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=  ~^p3(/i)COS 


(4.127) 
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1-  vdr 
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,(r,z) 


X.  nnz 
p4(n)Sin-— 


(4.128) 


1  d 


l-vdz 


(2-v)V2-|^ 

<7Z 


(r.z) 


Z=±L 


(4.129) 
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where  #,<„),  p2(«),  p3(„),  p4(n),  p5(?0f)  and  £>6«r0i)  are  given  by  Eqs.  (4.107),  (4.108), 

(4.115),  (4.116),  (4.119)  and  (4.120),  respectively.  Eqs.  (4.125)  through  (4.129)  specify 
six  boundary  conditions  for  £(r,Z),  however,  eight  boundary  conditions  are  required  for 
solving  Eq.  (4.121).  Hence,  the  two  needed  boundary  conditions  are  derived  from  the  fact 
that  the  shear  stress  does  not  exert  a  forcing  function  at  two  ends  of  the  coil.  From  Eq. 
(4.124),  we  may  write 

"  3"  (l-V)V2-—  §(,.«) 

1  -  v  dr  \  dz  ) 

The  substitution  of  §(r.«)  =  §,(r.z)  +  £2  (/•,*)  into  equation  (4.121)  yields  two  partial 
differential  equations  for  £,(r,Z)  and  £2(r,Z). 

V4^(r,z)  =  0  (4.131) 

V4£2(r,z)  =  0  (4.132) 

The  nonhomogeneous  boundary  conditions  for  ^r,Z)  do  not  fall  into  the  four  traditional 
(Cauchy,  Dirichlet,  Neumann  and  Robin)  categories.  Thus,  the  method  of  solving  the 
differential  equation  for  £(r,j>  is  somehow  different  from  the  usual  transient  method.  In  the 
usual  transient  method,  the  radial  nonhomogeneous  with  the  axial  homogeneous  boundary 
conditions  and  the  axial  nonhomogeneous  with  the  radial  homogeneous  boundary 
conditions  are  assigned  to  ^(r,Z)  and  £2(r,Z),  respectively.  Thus,  ^,(r.Z)  and  %2(r,Z)  each 
will  have  four  nonhomogeneous  boundary  conditions  in  one  direction  and  four 
homogeneous  boundary  conditions  in  the  other  direction.  Then  from  homogeneous 
boundary  conditions  and  using  one  of  the  conventional  methods,  e.g.,  eigenfunction 
expansion,  the  partial  differential  equation  is  reduced  to  an  ordinary  differential  equation. 
The  resulting  arbitrary  constants  from  the  solution  to  this  ordinary  differential  equation  are 
then  calculated  from  the  nonhomogeneous  boundary  conditions.  Hence,  solutions  for 


=  0. 


■=+/. 


(4.130) 
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£,(r,z)  and  £2(r,z)  are  obtained  independently  and  the  final  solution  for  £<m)  is  determined 
from  the  superposition  of  these  two  solutions. 

However,  In  the  present  formulation,  equations  for  the  boundary  conditions  are 
very  complicated  and  as  a  result  it  is  not  possible  to  reduce  the  partial  differential  equations 
to  ordinary  differential  equations  by  using  the  homogeneous  form  of  these  boundary 
conditions.  Therefore,  by  assuming  homogeneous  Dirichlet  boundary  conditions  in  the 
radial  and  axial  directions  for  £,(r,z)  and  £2(r,Z),  respectively,  the  partial  differential 
equations  are  reduced  to  ordinary  differential  equations.  Solutions  to  these  ordinary 
differential  equations  are  obtained  without  applying  the  nonhomogeneous  boundary 
conditions  and  they  are  superposed  to  give  the  solution  for  §(r.«).  Thus,  the  solution  for 
£(r,*)  will  have  eight  arbitrary  constants,  which  can  be  evaluated  from  the  eight  original 

boundary  conditions,  Eqs.  (4.125)-(4.130). 

The  partial  differential  equation  for  £,<r,z)  is  assumed  to  have  homogeneous 
Dirichlet  boundary  conditions  in  r.  Thus,  by  using  the  finite  Hankel  transform  (of  order 
zero)  in  the  radial  direction,  the  partial  differential  equation  for  |,  (/■.*)  is  reduced  to  an 
ordinary  differential  equation.  Applying  the  finite  Hankel  transform  to  Eq.  (4.131)  yields 

%[^]  =  {-a+fj)  £«..«>- 0  (4.133) 

where  £,(f0i,z)  is  the  transform  of  the  £,(<-, z>. 

b 

6  <f«    =  J  r&  c  ^o<f«  •  r^dr  (4. 1 34) 

Eq.  (4.133)  is  a  linear  ordinary  differential  equation  in  z  with  constant  coefficients. 
Solution  to  this  differential  equation  may  be  expressed  by 

J  (c0, ,  z)  =  AeioiZ  +  Be'^z  +  Cze(ojZ  +  Dze~Uz  (4. 1 35) 
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where  A,  B,  C  and  D  are  arbitrary  constants.  By  taking  the  inverse  finite  Hankel  transform 
of  Eq.  (4.135),  the  solution  for  £,o-.*)  is  obtained. 


(4.136) 


+C,z^wZ  +  £).ze-fo'z]^0(foi,r) 


The  partial  differential  equation  for  £2<r,*)  has  homogeneous  Dirichlet  boundary 
conditions  in  z.  Thus,  the  finite  Fourier  sine  transform  in  axial  direction  is  used  to  reduce 
the  partial  differential  equation  for  £2(r.«)  to  an  ordinary  differential  equation.  Hence,  by 
applying  the  finite  Fourier  sine  transform,  we  may  write 


3S[V4&(r.z)]  = 


(  X 


2_2  \ 


d2     1  d  n2K 
+  — — - 


ydr2    rdr  L 


£2(r.n)  =  0 


(4.137) 


where  £2(r.»)  is  the  finite  Fourier  sine  transform  of  the  function  £2  (r.«). 


£2(r,;i)  =  J  £2(r,z)Sin  


(4.138) 


Eq.  (4.137)  is  a  linear  ordinary  differential  equation  in  r.  Solution  to  this  equation 
may  obtained  by  introducing  a  new  function  S(r,»). 


S(r,n)  = 


'<92    1  <?  «VV 

H  Z  —  £2(r,n) 


dr2    r  dr  L2 


(4.139) 


J 


Incorporating  Eq.  (4.139)  to  Eq.  (4.137),  yields  a  differential  equation  for  H(r.«) 


(  32 


d2     1  d  n2K 


2  mi  \ 


ydr2     rdr  L2 


E(r,n)  =  0 


(4.140) 
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Eq.  (4.140)  is  a  modified  Bessel  equation  of  order  zero.  Solution  to  this  equation  is  given 
by 


S(r,»)  =  AI0&r)+  BK0&r) 


(4.141) 


where  A  and  B  are  arbitrary  constants,  and  70(xr)  anc*  ^o(xr)  3X6  tne  modified  Bessel 
functions  of  the  first  and  second  kind  (order  zero),  respectively.  The  substitution  of  Eq. 
(4.141)  into  Eq.  (4.139)  gives  a  differential  equation  for  £2(r.n). 


( d2     Id  n2n 


2  ~\ 


ydr2     rdr  L2 


%2(r,n)  =  AIQ(~r)  +  BK0(~r) 


(4.142) 


The  homogenous  and  particular  solutions  to  Eq.  (4.142)  may  be  written  as 


£2  (r,  n)  =  C/0       r)  +  DKQ(™  r)  +  F/7,  (™  r)  +  FrAT,  ( ™  r) 


(4.143) 


where  C  and  D  are  arbitrary  constants,  and  E  and  F  are  constants  that  must  be 
determined  from  the  nonhomogeneous  part  of  Eq.  (4.142).  Substituting  Eq.  (4.143)  into 
Eq.(4.142)  results  in 


r  d*  +L4--^P\(ErIl^r)  +  FrKi («„)-  AI0^r)  + BK0^r).  (4.144) 


dr2    r  dr  Is 


Expanding  the  left  hand  side  of  Eq.  (4. 144),  leads  to 


dr2     r  dr  Is 


(  # 


d2     1  d  n2K 


dr2  +  rdr  L2 


1  r) 

A<-rr>  +  —  I\(-rr)  +  2  —  L(-rr)E 

r  dr 


+ 


(  J2 


<?z     1  (9  n2K 


2_2  \ 


v<?r2  '  r<9r  L2 


,  1      me  d 


(4.145) 


ff1(£r)+-tfI(£r)  +  2-^-A:i(£r)F 

r  dr 
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From  the  modified  Bessel  equation  we  may  write 


dr1 


1  d  n2n2^ 

+  r- 


dr 


) 


/nit  1    j  nit 

1(xr)  =  -j/1(-rr) 


and 


V 


dr2+rdr      L1  , 


nit  1    zr  nit 

Ai  (~T  r)  ~  ~ 2    l  (X  r)  • 


Introducing  these  expressions  into  Eq.  (4.145),  yields 


'  d2     1  d     n2K2  \  a  -  \ 


(4.146) 


7  1<T"r    ^  l<T"r) 


£  + 


r  dr 


By  using  identities  /;(£r)  = /o(£r)- and  K[^r)  =  K^^-^K^r)  Eq. 

—  r  — r 

L  L 

(4.146) reduces  to 


+  —.  r-  lErl^ri  +  FrK^r)) 


dr2     r  dr  L 


f 

<  L 


T    nit  1  r 


J] 


+ 


nK 


—  If    nit       .  Tr    nK  ■  v  M 


or 


«?r2 


r  or 


(4.147) 


-  Z  —  iiY,  (-j- r)  +  2  — A,  (— r). 

/  .  /  , 
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Incorporating  Eq.  (4.147)  into  Eq.  (4.144)  gives  expression  for  E  and  F . 


E  =  —A  (4.148) 
2nit 

A  La  * 

F  =  —B 
2nK 


Hence,  £2<r'n)  may  ^  written  as 


|  (r,n)  =  <:/  <-, +  + _^ A  rI^r)  +  —BrK.&r).  (4.149) 

1  L      Inn         L      2nK  L 


The  solution  to  £2(r,Z)  is  simply  the  inverse  finite  Fourier  transform  of  £2(r,n). 

=  Ss't^.")]  =  7X?2<r'">sin^  (4. 150) 


Introducing  Eq.  (4.149)  into  Eq.  (4.150),  gives  the  final  form  for  ?2(r.Z). 


1  °° 


(4.151) 


sin- 


nnz 


The  superposition  of  Eqs.  (4.136)  and  (4.151)  furnishes  the  solution  to  £(r,z>. 

fr.o-gi  /2CoV°2(y  (4.152) 

+C;z^^  +  Z),.ze-^]/s:o(foi.r, 


+ 

L 


sin  
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Function  |(r,z)  is  an  odd  function  in  z,  hence  §0-.*)  =  -£(r.-*).  Incorporating  this 
expression  into  Eq.  (4.152)  leads  to 

[,Vfo'z  +  fl^~fo"  +  qze^1  +  D.tzeUz  ]  (4.1 53) 
=         fo'z  +  B/0iZ  +  Cfze~foi*  -  DtZeCmt] 

or 

Ai=~Bi  (4.154) 

q  =  Dr 

Thus,  the  expression  for  £(r,Z)  reduces  to 

fr-o-^jL^'y    [Asinh(fe)  <4155> 

+  C,zcosh(C0iz)]K0(corr) 


sin  


where  sinh(£0/z)  and  cosh(£0(z)  are  hyperbolic  sine  and  cosine,  respectively. 

Eqs.  (4.125)  through  (4.130)  describe  the  boundary  conditions  for  £(r,Z).  In  order 

to  evaluate  the  boundary  conditions,  we  need  first  to  derive  the  following  expressions. 
Laplacian  of  ^(r,z>: 

V%r,*=x>±  5oVo2^)    qsmCoiZ)K  (4.156) 

1  V  T  A   T    «'     _l  D  V    "*    1   '  HKZ 
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Second  partial  derivatives  of  £o-,z)  with  respect  to  z: 


+  Cizcosh(£0iz)]K0(t;0i.r) 
+  2^7°(y  C,.sinh(C0,z)^0(C0,.o 

-X-^[Q7o<Xr)  +  A,*()<x'> 

<xr)J  ~zT 


Second  partial  derivatives  of  §<r,z)  with  respect  to  r: 


-*2X    V°(y  [Asinhfe) 

/=1  J0^~J0^ 

+  Cizcosh(C0iz)]K0((0i.r) 

b0iJ0  ^o," 


-^£y2Wo(y  [Asinhfe) 


1  a 


+  C;zcosh(C0iz)]— r-^o«« 
r  or 


2n;r 


+ — D  +  —  Bn rKx (~  r, 


2n;r 


sin  


nKz 


-If;[c./,<¥'»+AJf,<?o]sto 

n  =  l  ^    '  ^ 


(4.158) 


2L 


2,  [  A,  4>  <-  +  BnK0  (T  r)l  sin  — 


First  partial  derivatives  of  V2£(r,Z)  with  respect  to  z: 
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j-V'fr,.,,-*2!    2^f    C.coshfe)^,.,)  (4.159) 


n=l  ^ 

+  fi„*T0<£,)]cos^ 


First  partial  derivatives  of  V2£(r.Z)  with  respect  to  r: 


j-V^^jr  J&S^g    qsinhC^)^^..)  (4.160) 


+  B„A:i(1-r)]sin-— 


Third  partial  derivatives  of  £<r,Z)  with  respect  to  z. 


jr^t)  =  ^±  [Acoshfe) 

+  Qzsinh(C0iz)]Ko(i;Oi,r) 

+  *2tJaj2oilf  C,cosh(C0,z)/r0(c0,,r, 
,=i  Jo^-Jo^ 


n  =  l 


L  *  t  nn  L  *  „„  nnz 
2niz         L      Inn  1  L 


(4.161) 


First  partial  derivatives  of  £(r.Z)  with  respect  to  r  and  second  partial  derivative  with  respect 
to  z: 
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~  fa*  =  ^I.^V  [A  sinh(  ;0iz) 
drdz  ^  J0((0ia)-J0a0ib) 


+Cizcosh(£0jz)]— K^.r) 


+  *2±  5°'J°y  C,.sinh(C0,z)f 
,=i  Jo^-Jo^b)  dr 


r) 


sin  


(4.162) 


First  partial  derivatives  of  ^<r,Z)  with  respect  to  z  and  second  partial  derivative  with  respect 


to  r. 


~fa*  =  -^t&ly  [M  cosh(M 

dzdr  £f  J0a0ia)-J0(c0ib) 


-4  ,2 

'Pro ( 

+  C,.(Co(zsinh(C0/z)  +  coshfe))]/^., 

-2  72 


(4.163) 


r) 


,=1  70(f0ia)-70(f0i*) 


1  <? 


+  C;(Co,zsinh(Co,z)  +  cosh(C0lz))]-— 

Jr  or 


■<  n37T3  r 


cos- 


n=l     ^  r 


nTtz 


tlKZ 


Note  the  difference  between  the  Fourier  Bessel  kernel  K0((0..r)  with  two  arguments,  and  the 
modified  Bessel  function  of  the  second  kind  K0(™r)  with  one  argument.  Now,  using 
K0(t0..b)  =  K0icora)  =  0,  the  following  evaluations  are  performed. 
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Eq.  (4.159)  at  r  =  a  : 


TV2^> 
oz 


=X^[  Vo<x^)  +  ^o<T-a) 


cos- 


nKz 


(4.164) 


Eq.  (4.161)  at  r  =  a: 


dz 


«  7T 


(4.165) 


cos- 


Eq.  (4.159)  at  r  =  fc: 


^2 


_Y>  nit\ 


nitz 


,[AnI0^h)  +  BnK0^h)]cos 

r=b      n=l  L  L 


(4.166) 


Eq.  (4.161)  at  r  =  b  : 


n  K 


r=b  »=1 


(4.167) 


cos- 


From  the  equations 


2  •/0(C0,fc) 
fl»  70(f0i.a) 


and 


r  =  /> 


_2_ 


at  r  =  a,  Eq.  (4.160)  becomes 
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+  2.7^[A«/i(-a)  +  B«^(xu)Jsm-T' 


(4.168) 


and  at  r  =  b  it  is 


or 


,<0 


•C,sinh(C0,z) 


+  2,    [  V.  <-r fc)  +  B«*i  (x  fe)J sin 


(4.169) 


With  the  same  substitution  at  r  =  a  Eq.  (4. 162)  generates 


dr  dz  : 


■fa 


=  ~t  ffi^^T^A  sinh(Co;z)+C,zcosh(C0,.z)]  (4.170) 


«v&&toC(Sinh(;(Z) 


a~i  J, 


En  /( 
r4 

(i=l  ^ 


sin- 


tlKZ 


and  Eq.  (4.163)  gives 


dz  dr 


=  4i^^[CAcosh(C0,,) 

+  C,(Co,zsinh(C0,z)  +  cosh(C0,z))] 

F!=l  ^ 


+^AnaIl^a)  +  -^-BnaKl^a) 
2h;t         L      2«7T  L 


cos- 


+  2,T^  Vohra)  +  BA<xa)  cos— 


nnz 


(4.171) 
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Moreover,  exploiting  the  same  substitution,  at  r  =  b  Eq.  (4.162)  transfers  to 


dr  dz 


r  =  h 


=  -fl  /22C°'J°(y)    [Asinh(C0,z)+C,zcosh(C0,Z)]  (4-172) 

0  ,=,  J0((0ia)-J0(S0ib) 

n  =  l  ^ 


2/1 7T 


2njz 


sin- 


and  Eq.  (4.163)  becomes 


JLi! 


■£(r.*) 


= 4i  y22^y  m  (4. 173) 

+  C,(Co,zsinh(C0,z)  +  cosh(C0,z))] 

K  71      X  .    nn      t   r\  v  Ml 


L     2  U J   1*  L    A  uv.  nx  1 


v  »2;r2  1  r  A  ,  n*     .  A  v  nx  1 

+  2,T7T  A, W»  +  fl«^o(-r«Jcos— . 


Now,  letting  z  =  L  results  in  the  following  equations. 
Eq.  (4.159): 


dz 


r) 


Eq.  (4.160): 


(4.174) 


=         J^'Y  C,.cosh(C0,.L^0(c0, 
+I^[Vo(x')+^o(x'>](-l)" 

n  =  l  ^ 


=  *T //^ C,sinh(f0,X)|g0,{„„  (4.175) 
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Eq.  (4.161): 


dz 


WW  [ACOSh(C0,X) 


^  702(f0ia)-/02(f0jfe) 
3  3 

w  7T 


+C;Lsinh(Co,X)]^0(fo,, 


CfCOSh(^)|X)K0(f(,.r) 


n=l  ^ 


2n7T 


2«;r 


(-1)" 


(4.176) 


Eq.  (4.162): 


dr  dz' 


Z=L 


=  ^If2U(y  [Asinh(C0,X) 


+C,Lcosh(£0iL)]— tfo(c0,.r) 


+         /22C°'Jo  T    C,  smh(COiL)fKo«0„r) 

£7  7n  (£,«)-  7n(CJ»  Cfr 


(4.177) 


Applying  the  shear  stress  boundary  condition  at  z  =  L  and  by  using  Eqs.  (4.175) 
and  (4. 177),  we  may  write 

m  70         70(fw«  <?r 


a)  —  J 


-  7T2I  -2  [A  sinh(Co,L)+C,.Lcosh(Co,L)]f 

/=i  70  (?„,«) -7o(?0,fc)  c?r 


"70        70(caj*)  ofr 


or 


,  =  1  »/( 


„  C°'J°(f;;>     {2  vC,  sinh(C0,7)  (4. 179) 

+  Co,[A  sinh(Co,L)+C,7cosh(C0,L)]}|:  K^.r)  =  0. 
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The  summation  in  Eq.  (4.179)  is  zero  if 


2  vQ  sinh(£0,L)  +  Co,[A  sinh(C0,L)  +  C,Lcosh(£0,.L)]  =  0  (4. 1 80) 


or 


[2v+C0iLcoth(C0,L)]C,.  +  C0,A=0. 


(4.181) 


Hence,  C(  is  expressed  by 


q  =  ft)(fo;.)A,. 


(4.182) 


where  coth(£0,L)  is  hyperbolic  cotangent  and  coaoi)  is  defined  by  Eq.  (4.183). 


coa  )  =  —  

2v+C0iIcoth(£0,I) 


(4.183) 


Substituting  Eqs.  (4.164)  and  (4.171)  into  the  radial  stress  boundary  condition  at 
r  =  a  and  using  the  orthogonality  property  of  cosine  yields 


(  \\nit 


v  — 

I  2) 


2  2 


L 
1 W 


2  h(Ta^~ 


j  nn 

2L3  ,(T0) 


V  

^  2J 


2  2 

nn  fl   Tt        rr  n% 


B,, 


+ 


n3K3 


3  3 

n  K 


r    nit  . 

I0(-ra)  + 


nV  1  ,  nn 

— i  V-rfl> 


C„ 


r4  "0lT 


nVl 

L3  a 


V  Ml 


D„ 


*  V  2CoVq(^)>/o(^) 


+  G 


Co,  A  fcos— cosh(^,,z)d!z 
-z,  L 

z. 

£0(.  Jzcos-— sinh(^0(z)Jz 
-z,  L 

+  fcos^^-cosh(^0(z)^z 


(4.184) 


=  -£>,(„)(!-  v). 


The  integrals  in  Eq.  (4.183)  are  evaluated  from  Eqs.  (4.185)  and  (4.186). 
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f  WfonttCrfVfe  =     "  (       2  - n V)sinh(C0lX)  (4.185) 

+  (C02,L2  +  «V)C0,Lcosh(£0l.L)]  =  O,  <*..«> 


jcos^cosh(C0iz)*  =  ?"sm{otL)  =  Q 

-L      1  \£oiL  +n  7t  ) 


2(f0i.,n)  (4.186) 


Incorporating  Eqs.  (4.182),  (4.185),  (4.186)  into  Eq.  (4.184),  gives 


A  A.  A  y\  ■  ^ 

Xi,t»)Ai  +  ^12(">fi«  +  fcl3<»>Q  +  ^14  <»>A  +  2^  Al(f0,-">A-  =  #1<»>  (4. 187) 

i=l 


where  Xuw,  X12(«),  X13(«)  and  X14<«)  are  given  by  Eq.  (4.188) 


X,|(n)  — 
X|2(n)  = 

X,j(B)  = 
X14(n)  = 


-1 


(i-v)L 
-l 


(l-v) 

-i 


r  n 

v  — 

I  2j 


2  2 

/itt  "  ™       7  nl 

°(xfl)~Tzr  i<xfl) 


(1-V) 
-1 


(1-v) 


—  /o<X«>  +  -73--'l<T-«> 

^    IM       ,   «2^2   1  v  nil 
 ^4_A0(Xa)  +  -^3_~Al(Xa) 


(4.188) 


and  AjCf^.n)  is  defined  by  Eq.  (4.189). 


_J  *  2Co,-/o(C„.a)70(g0,ft) 

(i  -  v)  a  y0  (f0.fl)  -  y0  ((„.*) 


[CoA(f«rf.") 


(4.189) 


+ ®<c«<>(CoA««-«) + Q2<^-."))] 


Enforcing  other  radial  boundary  conditions,  Eqs.  (4.126),  (4.127)  and  (4.128),  leads  to 


X21(n)A„  +X22(n)Bn  +X23(»)C„  +%24(n)Dn  +  ^A2((m,n)Ai =  p2w  (4.190) 


i=l 


*  See  Appendix  F  for  details. 
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^31<B>Ai  +%n(n)Bn  +X33(n)Cn  +%iA(n)bn  +       Aj^.n*^  =  p3(n)  (4.191) 


,'=1 


X4,(«)y4„  +X42(n)Bn  +X43(n)C„  H-X^  («)£>„  +XA4(^or")A  =P4(«)  (4.192) 


where  X21(»),  X22(n),  X23(«),  X24(«)  and  A^c^.io  are  given  by  Eq.  (4.193)  and  (4.194), 


X2|  (n)  — 


-1 


X22(n)  — 
X23(n)  = 


(1-v) 

-1 


(1-v) 

-1 


v  — 

/ — 

2j 


X24(n)  — 


(1-v) 

-1 


jik    nK      n2K2  nK 


2  2 

nr  W  »    i  v  nit 


(1-v) 


.—/0(X6)+— -/,(_/,) 


»,3_3  „2_2  i 

-p-A0(-r6)  +  -p--^1(-rW 


(4.193) 


(l-V)b  J0(C0ia)-J0(;0ib)L 


(4.194) 


+       {Co  A        +  ^2  (f«  •  »>)]. 


and  X31(n>,  X32(«),  X33(n>,  X34(n>  and  A3«;0/,n)  are  expressed  by  Eq.  (4.195)  and  (4.196), 


A31(«)-- 
X32(«)  = 

X33(n)  =  ' 


tin    nK         1  n2;r2 

7,(-;-a)  +  —  _,  ClLi-i-a) 


L2  1  1      (1  -  v)  2£ 

nit      nK            1      n27T2  rur 
— ^-A,(-j-a)i  5-  <ZA.n(-=-a) 

L     1  (1-v)  2L3  0 


3  "'0lX 


(4.195) 


nV  1 


,(-T-") 


L4  (1-v)  1:1 


X34(n)  —  — 


»V  1 


V  (1-v) 
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(4.196) 


and  finally  A4({0<.»)»  X41(«),  X42(«),  X43(«),  and  X^u)  are  defined  by  Eqs.  (4.197)  and 
(4.198). 


.  1     K     2QJ0  {{..aj)     r.  - 

l- y0(?0i.fl)-y0(f0,fc) 


(4.197) 


X4)(n)  —  — 


X42(«)  —  — 


X43(n)-- 


X^*")  — 


1  n2K2 


nit  nK 
—z-l,  (-rb)  + 

L   1         (1-  v)  2L3 


1  nV 


(l  -  v)  IE 


nV  1 


L   (1-v)  1 


nV  1 

L4  (1-v) 


(4.198) 


The  definitions  for  QjCf^.n)  and  Q4(f0,,n)  are  given  by  Eqs.  (4.199)  and  (4.200), 
respectively. 

jZsin^cosh(£0,z)<fe  =     2L^~lJ   2[2nxL£0ismh(C0iL)  (4.199) 

-  [Qfi  +  nV)n7rcosh(£0,.L)]  =  Q^.n) 


fsin^sinh(£0lz)<fc  =  jw0r  r3Tsinh^o,L)  =  Q^.-)  (4.200) 


-L 
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By  applying  the  axial  stress  boundary  condition  and  using  Eqs.  (4.174)  and 
(4.176),  we  may  write 

2£o/[  (1  -  2  v)ft)(fo.)Cosh(C0,L)  -  Co,  cosh(£0,X)  -  C0,G)(fo,.)Lsinh(C0,L)]A,  (4.201) 


+Z(-i)'f 


A, 


+S(-i) 


n=l 


-S(-i) 


n=l 


(2-v)|r^0(fw.o/0<?r)dr+||/r2J8:0(f«.o/1(^^r 

a  a 

(2-  V^r^^.^^^r  +  ^JrX^.^^^r  |fin 

a 


=  -(l-v) 


£7 

Jr^0(f0i,r)/s:0(^r)rfr 


*  Co; 


1- 


'  1  V 


By  introducing 


A5(fo,)  =  2Co2,[(1  -  2  v)(O(fO(.)Cosh(C0,L) 

-  Co,  cosh(C0;X)  -  Co,«(?0,)^sinh(C0iL)] 


(4.202) 


and 


*52<C«.»>  =  H)  "75" 


(2-v)J^o(fo,.,r)/o(^^r  +  ||Jr2/i:o(,0,.r)/1^^r 


(4.203) 


(2-v)Jr^0^,^o(x^  +  ffJ^o<fo,-^1(x^ 


^53^o,'»)  =  (-1)  -74- 


^54(fo,-'"  =  (-1) 


(7 

Jr/i:0(f0,,r)/0(^r)^r 

.a 
"fc 

jrKoa0i,r)Ko(^-r)dr 
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r5(fw.fli.)=-(i-v) 


1- 


Eq.  (4.200)  is  reduced  to 


(4.204) 


+ x53(f0,,«)C„  +  x54(f0(,n)Dn  =  r5(?0i.?li) 


(4.205) 


Eqs.  (4.187),  (4.190),  (4.191),  (4.192)  and  (4.205)  represent  a  system  of 

A  A  A  A 

equations  where  unknowns  are  An ,  Bn,  Cn,  Dn  and  A, .  As  it  is  practiced  in  numerical 
simulations,  the  infinite  series  is  changed  to  a  finite  summations  with  an  acceptable 
truncation  error.  Hence,  the  upper  limit  for  i  and  n  are  changed  to  M  and  N,  respectively. 
Expanding  these  finite  summations  would  result  in  a  system  of  equations  with  47V  +  M 

A  A  'A   '.         A  A  A  a  A. 

unknowns  and  equations.  Unknowns  are  A,---^,  Bl  -  Bn,  Cx  -  Cn,  Dx  -  Dn  and 
A,  -  •  •  j4n.  Eq.  (4.204)  gives  M  equations  by  letting  i  vary  from  1  to  M.  Moreover,  allowing 
n  to  advance  from  1  to  iV  in  Eqs.  (4.187),  (4.190),  (4.191)  and  (4.192),  results  in  4iV 
equations.  The  coefficient  matrix  for  this  system  of  equations  is  shown  in  Fig.  4.2.  By 
solving  this  system  of  equations,  the  arbitrary  constants  for  the  complementary  solution  of 
the  axial  stress  function  are  obtained.  The  combination  of  the  complementary  and  the 
Green's  function  solutions  for  the  axial  stress  function  yields  a  solution  that  satisfies  both 
the  boundary  conditions  and  the  axial  body  force.  This  solution  together  with  the  radial 
Green's  function  determine  the  disturbance  of  stress  and  strain  in  the  given  coil. 
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CHAPTER  5 
RESULTS  AND  DISCUSSIONS 

Introduction 

Results  of  the  generalized  plane  strain  analysis  and  Green's  functions  solution  are 
presented.  A  comparison  of  the  analytical  results  with  the  finite  element  results  is 
performed  to  check  the  validity  of  the  solutions.  The  finite  element  results  are  obtained* 
through  the  ABAQUS30  finite  element  solver.  Linked  with  the  finite  element  solver  is  an 
in-house  FORTRAN  subroutine31  which  determines  the  non-uniform  body  forces  created 
by  the  magnetic  fields.  These  forces  are  calculated  at  each  gaussian  integration  point  and 
are  used  as  an  input  to  the  finite  element  program. 

Results  of  Generalized  Plane  Strain  Analysis 

An  example  magnet  is  generated  to  examine  the  results  of  a  generalized  plane  strain 
analysis.  The  example  design  is  for  a  15  Tesla  magnet  with  two  coils  each  with  two  layers 
(Fig.  5.1).  The  parameters  for  each  of  the  coils  are  given  in  Table  5.1,  where  for  each 
layer  the  inside  radius  is  a\,  the  outside  radius  is  a2,  the  half  length  is  L,  and  the  current 
density  is  J.  The  outer  radius  of  the  reinforcement,  where  applicable,  is  given  by  a$. 


Table  5.1  Parameters  for  15  T  example  magnet  with  two  coils,  each  with  two  layers. 


Coil 

a\  (mm) 

CJ2  (mm) 

as  (mm) 

L  (mm) 

J  (A/mm2) 

la 

100 

126 

250 

102.36 

lb 

126 

148 

154 

250 

161.37 

2a 

157 

170 

250 

155.55 

2b 

170 

210 

250 

149.54 

Finite  element  analysis  was  done  by  Iain  Dixon  at  National  High  Magnetic  Field 
Laboratory. 
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The  example  magnet  has  two  combined  layers  containing  Nb3Sn  conductor  in  the 
inner  coil,  and  two  connected  layers  containing  NbTi  conductor  in  the  outer  coil.  The 
second  layer  of  the  inner  coil  contains  a  reinforcement  on  the  outer  diameter,  as  indicated  in 
the  Fig.  5.1. 


Figure  5.1  15  Tesla  example  magnet  with  two  coils,  each  with  two  layers. 

The  difference  in  conductors  is  reflected  in  the  coil  material  properties  listed  in 
Table  5.2.  It  is  noted  that  the  properties  in  the  tangential  direction  are  dominated  by  the 
properties  of  the  conductor.  The  properties  transverse  to  the  conductor,  in  the  r  and  z 
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directions,  are  influenced  strongly  by  the  properties  of  the  insulation  between  the 
conductors. 


Table  5.2  Materials  properties  of  coils  and  reinforcement. 


Nb3Sn  Coils 

NbTi  Coils 

Reinforcement 

kg  (GPa) 

66  500 

1 90  000 

br  (GPa) 

41  400 

45  400 

1 90  000 

bz  (OPa) 

44  500 

49  400 

1 90  000 

14  700 

iTi  /  \J\J 

16  000 

1  U.vV/v 

7^  000 

G.  (GPa) 

15.700 

17.500 

73.000 

G„  (GPa) 

11.800 

12.500 

73.000 

0.341 

0.341 

0.305 

0.318 

0.319 

0.305 

v„ 

0.204 

0.182 

0.305 

aeAT 

-0.0028 

-0.0028 

-0.0029 

arAT 

-0.0048 

-0.0048 

-0.0029 

azAT 

-0.0048 

-0.0048 

-0.0029 

Magnetic  Stress  and  Strain 

The  following  is  a  comparison  between  results  generated  by  the  generalized  plane 
strain  and  the  finite  element  methods.  Stresses  and  strains  are  calculated  for  the  example 
magnet.  Moreover,  plane  stress  and  plane  strain  analyses  are  performed  for  the  example 
magnet  and  compared  along  with  the  generalized  plane  strain  analysis. 

In  reference  to  Fig.  5.2,  the  tangential  stress  is  graphed  from  the  results  of  the  finite 
element  and  the  generalized  plane  strain  analyses.  Small  differences  are  displayed  with  the 
greatest  discrepancy  within  the  reinforcement.  This  may  be  caused  by  a  stiffer  solution 
resulting  from  relatively  fewer  elements  being  used  in  the  narrow  reinforcement  region. 

Fig.  5.3  shows  the  characteristics  of  radial  stress  through  the  magnet.  The 
difference  between  generalized  plane  strain  and  finite  element  for  radial  stress  is  even 
smaller.  Fig.  5.4  shows  the  characteristics  of  axial  stress.  Since  in  the  finite  element 
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methods  the  axial  strain  is  not  being  constrained  to  a  constant,  the  results  of  the  finite 
element  solution  is  improved  over  the  generalized  plane  strain  analysis. 

Figs.  5.5  and  5.6  are  plots  of  the  tangential  and  radial  strains.  Both  graphs  show 
an  extremely  close  relationship  between  the  two  methods.  Fig.  5.7  shows  the  variation  of 
axial  strain  through  each  coil  for  finite  element  along  with  the  constant  value  of  the  axial 
strain  for  generalized  plane  strain  analysis.  The  maximum  discrepancy  for  axial  strain 
between  these  two  methods  is  nine  percent. 

The  theory  of  generalized  plane  strain  is  compared  to  other  analytical  methods  as 
well  as  finite  elements.  Using  the  example  magnet,  the  results  of  the  generalized  plane 
strain  analysis  are  studied  along  with  those  determined  by  plane  stress  and  plane  strain 
assumptions. 

The  stress  in  the  tangential  and  radial  directions  for  the  three  analytical  solutions  are 
plotted  in  Figs.  5.8  and  5.9.  Results  match  quite  well  and  are  indistinguishable  through  the 
superconducting  layers.  The  reinforcement  possesses  the  largest  difference  between  the 
generalized  plane  strain  and  the  other  theories. 

The  greatest  improvement  in  the  approximation  of  stress  is  shown  in  Fig.  5.10, 
where  the  axial  stress  for  plane  strain  and  generalized  plane  strain  is  plotted.  Note  that  for 
the  plane  stress  condition  the  axial  stress  component  is  zero.  With  a  more  general 
approximation  of  axial  strain,  the  generalized  plane  strain  theory  is  able  to  predict  axial 
stress  to  a  higher  degree  of  accuracy.  Referring  to  Fig.  5.4,  it  is  easily  seen  that  the 
generalized  plane  strain  matches  the  finite  element  more  closely  than  the  plane  strain 
method.  In  fact  the  plane  strain  result  predicts  the  axial  stress  in  the  tensile  regime, 
whereas  the  generalized  plane  strain  and  finite  element  analyses  show  the  axial  stress  in 
compression.  From  the  direction  of  the  radial  and  axial  fields  produced  from  the  magnet,  it 
may  be  easily  found  that  the  axial  stress  is  compressive  in  nature. 

The  tangential  and  radial  strains  in  Figs.  5.1 1  and  5.12  respectively  show  the  same 
type  of  relationships.  The  strains  start  to  deviate  significantly  with  the  plane  stress  results 
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being  bound  by  the  other  theories.  The  axial  strain  for  plane  stress  and  generalized  plane 
strain  is  plotted  in  Fig.  5.13.  The  plane  stress  predicts  much  lower  values  for  axial  strain. 
Note  that  for  the  plane  strain  condition  the  axial  strain  component  is  zero. 

Thermal  Stress  and  Strain 

In  the  previous  comparisons  magnetic  forces  are  considered  as  the  applied  loads. 
Since  the  generalized  plane  strain  assumption  may  be  applied  to  the  thermal  as  well  as 
magnetic  loads,  this  case  is  also  studied  to  observe  the  effects  of  cryogenic  temperatures  on 
superconducting  magnets.  The  following  studies  the  same  15  Tesla  magnet  by  comparing 
the  stress  and  strain  results  determined  by  the  generalized  plane  strain  analysis  with  the 
finite  element  methods.  The  thermal  loads,  that  are  applied  in  terms  of  ccAT,  are  defined  in 
Table  5.2. 

The  comparisons  that  follow  are  similar  in  structure  to  the  magnetic  load 
comparisons.  The  differences  between  the  finite  element  and  generalized  plane  strain 
methods  are  even  less  because  of  the  near  constant  behavior  of  the  axial  strain.  In  Figs. 
5.14  -  5.18,  the  tangential,  radial,  and  axial  stresses  and  tangential  and  radial  strains  show 
excellent  agreement.  The  slight  differences  that  do  exist,  occur  more  frequently  in  the  first 
coil  than  the  second.  The  reason  for  this  may  be  seen  in  Fig.  5.19  where  there  is  a  larger 
variation  of  axial  strain  through  the  first  coil  than  the  second,  which  is  near  zero. 

Results  of  Green's  Functions  Solution 

The  Green's  functions  solution  is  applied  to  a  23  T  isotropic  coil  and  the  results  are 
compared  to  the  results  of  finite  element  methods.  The  parameters  for  this  coil  is  given  in 
Table  5.3,  where  the  inside  radius  is  a,  the  outside  radius  is  b,  the  length  of  the  coil  is  2L, 
the  current  density  is  J,  and  the  elastic  modulus  and  the  Poison's  ratio  for  the  coil  are  E  and 
v,  respectively. 
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Figs.  5.20  and  5.21  show  the  tangential  and  radial  stresses  through  the  coil  along 
the  radius  at  three  different  axial  position  (midplane,  half  between  the  midplane  and  the  top, 
and  the  top  of  the  coil).  Little  differences  are  displayed  between  results  of  the  Green's 
functions  solution  and  finite  element  methods  with  the  greatest  discrepancy  at  the  top  of  the 
coil.  Fig.  5.22  shows  the  characteristics  of  the  axial  stress  through  the  coil  along  the  radius 
at  midplane  and  z  -  L/2.  Note  that  due  to  traction  free  boundary  conditions,  axial  stress  at 


Table  5.3  Parameters  for  the  23  T  isotropic  coil. 


Name 

Symbol 

Value 

Unit 

Inner  radius 

a 

100.00 

mm 

Outer  radius 

b 

136.50 

mm 

Half  length 

L 

28.00 

mm 

Elastic  modulus 

E 

111.00 

GPa 

Poison's  ratio 

V 

0.30 

Current  density 

J 

530.10 

A/mm2 

the  top  of  the  coil  is  zero.  The  Green's  functions  solution  shows  a  good  correlation  with 
the  finite  element  for  axial  stress.  Fig.  5.23  shows  the  shear  stress  at  z  =  L/2.  The 
maximum  difference  between  the  results  of  finite  element  methods  and  Greens  functions 
solution  for  the  shear  stress  is  about  6%.  Shear  stress  is  zero  at  z  =  L  due  to  the  boundary 
conditions,  and  it  is  zero  at  midplane  because  of  symmetry.  Figs.  5.24  -  5.26  are  plots  of 
the  normal  strains.  These  graphs  show  the  same  close  relationship  between  the  Green's 
functions  solution  and  finite  element  analysis. 
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CHAPTER  6 
SUMMARY  AND  CONCLUSIONS 


Analytical  closed  form  solutions  for  the  distribution  of  stress  and  strain  have  been 
developed  in  high  field  solenoid  magnets,  including  superconducting  magnets.  These 
solutions  are  presented  in  the  form  of  generalized  plane  strain  analysis  and  the  Green's 
functions  solution.  The  generalized  plane  strain  analysis  provided  an  analytical  formulation 
for  the  magnetic  and  thermal  state  of  stress  at  the  midplane  of  a  coil,  where  shear  stress  is 
negligible.  This  formulation  considers  the  effect  of  the  axial  body  force  between  the 
midplane  and  the  end  of  the  coil.  The  three-dimensional  stress  analysis  in  terms  of  Green's 
functions  provided  solutions  for  components  of  stress  (including  the  shear  stress) 
throughout  the  coil.  The  Green's  functions  method  permits  the  development  of  a  solution 
irrespective  of  the  type  of  the  field  or  its  distribution  within  the  coil.  Thus,  the  Green's 
functions  may  be  used  for  both  thermal  and  magnetic  stresses.  The  problem  was 
formulated  in  terms  of  stress  functions.  Green's  functions  were  derived  by  exploiting 
finite  Hankel  and  finite  Fourier  transforms.  Boundary  conditions  were  satisfied  by 
introducing  a  complementary  solution  for  the  axial  stress  function.  The  radial  Green's 
function  with  the  superposition  of  the  complementary  and  the  axial  Green's  function 
provided  a  comprehensive  analytical  solution  for  the  stresses. 

A  comparison  of  the  analytical  results  with  the  results  of  finite  element  methods  was 
presented  to  confirm  the  validity  of  the  solutions.  The  generalized  plane  strain  showed  a 
good  agreement  with  the  finite  element  methods  for  the  radial  and  tangential  stresses  and 
strains.  The  resulting  axial  stress  and  strain  from  the  generalized  plane  strain  analysis  were 
in  an  acceptable  range  (less  than  9%  error)  compared  to  the  finite  element  methods.  The 
generalized  plane  strain  analysis  was  also  compared  to  other  plane  theories.  The 
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generalized  plane  strain  improved  the  prediction  of  axial  stress  and  strain  significantly 
compared  to  plane  stress  and  plane  strain  theories.  The  resulting  stress  and  strain  from  the 
Green's  functions  solution  showed  a  good  correlation  with  the  finite  element  methods. 

The  Green's  functions  solution  provided,  for  the  first  time,  a  complete  analytical 
stress  solution  for  an  isotropic  coil.  This  solution  should  be  used  as  a  foundation  for  the 
stress  analysis  of  multilayer  magnets.  The  future  work  should  also  extend  this  solution  for 
an  orthotropic  coil.  The  use  of  the  Green's  functions  solution  is  not  restricted  to  magnets. 
It  can  also  be  used  for  other  elasticity  problems  with  physical  or  fictitious  body  forces. 


APPENDIX  A 

EQUATIONS  FOR  THE  CONSTANTS  OF  THE  THERMAL  STRESS  ANALYSIS 


In  this  appendix  the  axial  strain  and  arbitrary  constants  of  the  thermal  stress 
formulation,  based  upon  generalized  plane  strain  analysis,  are  computed  for  the  different 
kinds  of  boundary  conditions. 

The  radial  displacement  and  the  radial  and  axial  stresses  are  expressed  by  Eqs. 
(A.l),  (A.2)  and  (A.3),  respectively. 

u  =  D/  +  D2rk  +  Axe'z°'r  +  A2r  (A.l) 

Or  =  (C12  +  kC21)D/-x  +  (C12  - kC22)D2r-k-]  (A.2) 
+(C12  +  C22)A\£'"'  +  (C12  +  C22)A2  +  C23e'°' 
-Cx2aeAT  -  C22arAT  -  C23azAT 

oz={Cn  +  kC2,)D/-x+(C„-kC2,)D2r-k-x  (A.3) 

+(Cl3  +  Crh)A\e7  +  (^13  +  ^23)^2  +  C33ez 


,((« 
'-z 


-Ci3agAT  -  C23arAT  -  C}3azAT 


For  a  single  layer  coil  with  constant  current  density  and  uniform  material  properties, 
the  boundary  conditions  are 

ar  =  0  at  r  =  a,  (A.4) 


<7r  =  0  at  r  =  a2  (A  5) 


"2 


j27trazdr  =  0.  (A.6) 
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Evaluating  equation  (A.2)  at  the  inner  radius  (a\)  results  in 


a11D1+a12D2+fl13ef'=^1 


where  the  constants  are  given  by 

au  ={Cn  +  kC22)a]k- 

an=  (^i2  —  ^nrh 
aB  =  (^12 +  ^22)  A  +  Q3 

b\  =  ~{c\2  +  Q2M2  +  CnaeAT  +  C22arAT  +  C23azAT. 


-k-\ 

1 


Evaluating  Eq.  (A.2)  at  the  outer  radius  («2)  results  in 

a21D,  +  a22D2  +  a23e'°'  =  b2 

where  the  constants  are  given  by 

a22  =  (Q2  -  kC22)a 


k-l 

2 


a23  =  «13 

b2  =bv 


Integrating  Eq.  (A.3)  in  Eq.  (A.6)  results  in 


a3A+anD2+a33eT  =b3 


where  the  constants  are  given  by 
%  ={c»  +  kC2i) 


k  +  \ 


a      -  a 

a32  =  (^13  ~  ^C23)    2  +  j 

°33  =  [(^13  +  Q3M1  +  C33]    2    2  ~ 

^3  =  -[(C13  +  Cq)^  -  Cna6AT-  C23arAT-  C33azAj] 


a2  - 
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The  set  of  linear  Eqs.  (A.7),  (A.9),  and  (A.l  1)  may  be  written  in  matrix  form 


«11 

an 

a13 

A 

V 

a2] 

a22 

D2 

a32 

a33_ 

.£z_ 

A. 

(A.13) 


which  is  solved  in  the  standard  way  as  given  by  Eq.  (A.  14). 


A 

au 

°12 

A 

a22 

a 

A. 

_a3l 

an 

a 

3. 


(A.  14) 


For  a  single  layer  coil  with  a  reinforcement,  the  boundary  conditions  are  applied  to 
the  coil  and  the  reinforcement. 


cr*1'  =0  at  r  =  a, 


(A.  15) 


oj"  =  <7<2)  at  r  =  a2 


(A.  16) 


=  w<2)  at  r  =  a2 


(A.  17) 


<7<2)  =  0  at  r  =  a, 


(A.18) 


The  axial  force  equilibrium  applies  to  the  coil  and  reinforcement. 


"j27urazdr  =  0 


(A.  19) 


Evaluating  Eq.  (A.2)  at  the  inner  radius  (a\)  results  in 


auD,+anD2+ax5e'^  =b, 


(A.20) 


where  the  constants  are  given  by 


-k-\ 

*2 
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au  ={Cn+kC22)alk~l  (A.21) 

an  =  (^12  ~~  ^22)^1 

^5  =  (^12  +  +  Q3 

^  =  ~(ci2  +  Q2  K  +  CnagAT  +  C22arAT  +  C23azAT. 
Evaluating  Eq.  (A.2)  at  the  interface  between  the  coil  and  the  reinforcement  (#2)  results  in 

a2\ D\  +  «22  A  +  a23D'l  +  a24  A'  +  a25£r  =  ^2  ( A.22) 

where  the  constants  are  given  by 

«2i  ={Cn+kC22)a2k-x  (A.23) 
fljj  =  (C12  -  kC22)a2 

^23  =  —  (^12      ^  ^22)^2 

a24  =  —  (C12  —  &  ^22)2^2 

a25  =  [(^12  +  Q2)A  +  Q3]  ~  [(^12  +  +  Q3] 

^2  =  [(CT2  +  Q'2  K  -  C>X     -  C22a'rAT  -  C2\a[AT]  - 
[(C12  +  C22  )A2  -  CnaeAT  -  C22arAT  -  CnazAT\ 

and  the  unprimed  and  primed  quantities  refer  to  the  coil  and  the  reinforcement,  respectively. 
Evaluating  Eq.  (A.  1)  for  the  displacement  at  the  interface  results  in 

<HA  +  a32D2  +  «33D1  +  «34  D2  +  <*3&*  =  bi  (A.24) 

where  the  constants  are  given  by 

a31  =  a2  a32  =  a2k  (A.25) 


k'-\ 

*2 


a33  =  -a2 


-k' 


a35  =  Axa2  -  A[a2         &3  =  A^a2  -  A2a2 


Ill 


Evaluating  Eq.  (A.2)  at  the  outside  radius  of  the  reinforcement  (03)  results  in 

«43 D;  +  ajy2  +  a45ef  =  b4  (A  26) 


where  the  constants  are  given  by 

"«={C[2+k'C21  )o3*'-1  (A.27) 

fl44  =  (^12  —  k  Q2)a3 
fl45  =  (Q2      Q2  ) A  Q3 

^  =  -(c;2  +  c2'2  ^  +  C,'2a;A7 + C2'2ar'Ar + c;3cc'zat. 


Integrating  Eq.  (A.3)  through  coil  and  reinforcement  in  Eq.  (A.  19)  results  in 

«S1 A  +  a52  D2  +  °53  A'  +  *54  A'  +  ^S^"'  =  ^5  (A.28) 


where  the  constants  are  given  by 

k  +  l  t+I 

"w-fa+AQ,)^    ^  (A.29) 


U52  ~  (^13  ^23) 


a54=(q,3-^) 


a2-^'-ar'+' 

fc'  +  l 

„  -k'+i  „ 
a3      —  a2 


-k'  +  l 


2^2  _ 2    „ 2 


«55  =  [(cl3  +  c23)a  +  c33]^l + [(q3  +  c2'3)a,  +  c3'3]^-^ 

fc5  =  -[(CI3  +  C23)A2  -  C13aeAr-  C23arAr-  C33azAr]- 


ta2  -a, 


2 

-  [(CJ  +  C2'3 -  q3a'eAT  -  C'2,a'AT  -  C3'3az'Ar] 


^3-^2 


The  set  of  linear  Eqs.  (A.20),  (A.22),  (A.24),  (A.26),  and  (A.28)  may  be  written  in  the 
matrix  form 
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«I1 

aI2 

0 

0 

"A" 

V 

a22 

a23 

a24 

a25 

A 

&2 

a33 

fl34 

fl35 

A 

&3 

0 

0 

«43 

fl44 

«45 

A 

*>4 

«51 

«52 

«53 

«54 

«55. 

A. 

A. 

(A.30) 


which  is  solved  as  given  by  Eq.  (A.31). 


"A" 

«12 

0 

0 

a15" 

-1 

V 

A 

a22 

«23 

a24 

«25 

*2 

A 

«31 

«32 

«33 

fl34 

*3 

A 

0 

0 

«43 

«44 

*4 

A. 

.a51 

«52 

a53 

«54 

fl55. 

A. 

(A.31) 


For  a  general  multilayer  coil  in  layers)  the  boundary  conditions  are  applied  to  each 


layer. 


ct*0  =0  at  r  =  a, 


^')  =  ^+l)  at  r  =  ai+l;j  =  ln-\ 


at  r  =  aj+];j  =  \,n-\ 


a\n)  =  0  at  r  =  a 


H  +  I 


(A.32) 
(A.33) 
(A.34) 

(A.35) 


The  axial  force  equilibrium  also  applies  to  all  layers. 

$2KrGzdr  =  0 


(A.36) 


The  derivation  of  the  equations  for  the  constants  is  the  same  as  the  previous  cases. 
The  matrix  of  coefficients  for  the  system  of  2n+l  equations  and  unknowns  is  given  by  Eq. 
(A.37). 
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All 

«12 

0 

0 

0 

0 

o  ... 

a\  2n+l 

<h\ 

«22 

a23 

«24 

0 

0 

o  ... 

a2  2n+l 

«31 

«32 

«33 

«34 

0 

0 

o  ... 

a3  2n+l 

0 

0 

a43 

«44 

«45 

a46 

0  ••• 

a4  2n  +  l 

0 

0 

a53 

fl54 

a55 

°56 

o  ... 

fl5  2n  +  l 

0 


*2n+l  1 


0  a 


In  2n-l 


2n+l  2 


2n  2n 


2n+l  2n 


2n  2n+l 


2n+l  2n+l 


(A.37) 


The  solution  to  this  system  of  equations  is  obtained  by  matrix  inversion. 


APPENDIX  B 

PARTIAL  DIFFERENTIAL  EQUATIONS  FOR  THE  STRESS  FUNCTIONS 


In  this  appendix  a  detailed  derivation  of  the  partial  differential  equations  for  the 
stress  functions  are  given. 

For  an  axisymmetric  distribution  of  magnetic  body  forces  X(r,z)  acting  in  the 
interior  of  an  elastic  isotropic  coil,  the  equilibrium  equation  in  the  vector  form  is  given  by 
Eq.  (B.l). 

V.C  +  X  =  0  (B.l) 

The  constitutive  equation  relating  the  stress  tensor  to  the  strain  tensor  is  expressed  by  Eq. 
(B.2). 

a  =  Afr(£)I  +  2jU£  (B.2) 
The  displacement  vector  U(r,Z)  is  related  to  the  strain  tensor  as  given  by  Eq.  (B.3). 

e  =  |[Vu  +  VTu]  (B.3) 

Incorporating  Eq.  (B.3)  into  Eq.  (B.2)  yields  an  expression  for  the  stress  tensor  in 
terms  of  displacement  vector. 

a  =  A(V.u)I  +  ju[Vu  +  VTu]  (B.4) 

Introducing  Eq.  (B.4)  into  Eq.  (B.l)  gives 

AV.[(V.u)I]  +  /iV.[Vu  +  VTu]  +  X  =  0.  (B.5) 

From  continuum  mechanics,32  we  may  write 
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V.[(V.u)I]  =  V(V.u)  (B-6) 
V.(Vu)  =  V(V.u) 
V.(VTu)  =  V2u. 

Incorporating  the  expressions  form  Eq.  (B.6)  into  Eq.  (B.5)  gives 

(A+/z)V(V.u)  +  /xV2u  +  X  =  0.  (B.7) 

From  the  Helmholtz  theorem,  any  vector  satisfying  Eq.  (B.7)  may  be  resolved  into 
the  sum  of  a  gradient  and  a  curl 

u  =  V0  +  VxA  (B.8) 

where  (fHr.z)  is  a  scalar  potential  and  A(r,2)  is  a  vector  potential  such  that  V.A  =  0. 
Incorporating  the  displacement  vector  from  Eq.  (B.8)  into  Eq.  (B.7)  gives 

(A  +       V.(Vtf>  +  V  x  A)]  +  juV2  (V0  +  VxA)  +  X  =  0.  (B.9) 

Nevertheless,32  the  divergence  of  the  curl  of  any  vector  is  zero  V.(VxA)  =  0,  the 
divergence  of  the  gradient  of  a  scalar  is  the  Laplacian  of  a  scalar  V.(V0)  =  V2</>,  the 

Laplacian  of  the  gradient  of  a  scalar  is  the  gradient  of  the  Laplacian  of  a  scalar 
V2(V.0)  =  V.(V2</>)  and  finally  the  Laplacian  of  the  curl  of  a  vector  is  the  curl  of  the 

Laplacian  of  a  vector  V2(Vx  A)  =  Vx(V2A).  By  using  these  identities  Eq.  (B.9)  is 
reduced  to 

(A  +  2//)  V(  V20)  +  /iV  x  (V2A)  +  X  =  0 .  (B.  10) 

The  independent  potential  functions  0  and  A  may  be  expressed  as 

A  =  aVxvP  (B.ll) 

where  a  and  (5  are  arbitrary  constants  and  the  vector  *F  is  the  stress  function  vector. 
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Introducing  Eq.  (B.l  1)  into  Eq.  (B.10)  gives 


P(X  +  2ju)V[v2(V.vF)]  +  a/iV  x  [V2(V  x  ¥)]  +  X  =  0  (B.12) 


or 


P(X  +  2ju)V[V.(V2XF)]  +  a/iV  x  [V  x  (V2VF)]  +  X  =  0 .  (B.  13) 


Substituting  V[V.(V2VF)]  in  Eq.  (B.13)  with  V2(V2VF)  +  Vx[Vx(V2XF)]  yields  Eq. 

(B.14). 

P(X  +  2/x)V4XF  +  [fi(X  +  2/x)  +  jia]V  x  [V  x  (V2XF)]  +  X  =  0  (B.  14) 

The  second  term  of  Eq.  (B.  13)can  be  eliminated  by  making  [fi{X  +  2/j,)  +  fia]  to  be 


zero.  The  first  term  of  Eq.  (B.13)  can  be  simplified  by  choosing  /?  as 


,  which 


results  in  a  equal  to  — - .  Replacing  a  and  B  respectively  with  — -  and  — - —  in  the 

fi  fi        A  +  2\l 

Eq.  (B.14)  yields  a  partial  differential  equation  for  the  vector 


V4VF  +  X  =  0 


(B.15) 


Rewriting  Eq.  (B.15)  in  component  form  yields 


V>-±T«F-±^.-±|  V2- 


-2  '    '    r4  dO2  r2 


r2J 


dd 


(B.16) 


V2-4rl 


4  d2Ya  4f 


dd2  V 


v 


n 

r2J 


V2- 


dd  6 


V4*Fz+Xz=0. 


Eq.  (B.16)  defines  the  three  partial  differential  equations  of  stress  functions. 


APPENDIX  C 

ORTHOGONALITY  PROPERTY  OF  THE  FOURIER  BESSEL  KERNEL 


In  this  appendix  the  orthogonality  property  of  the  Fourier  Bessel  kernel        r)  is 
proved  and  the  expression  for  the  lA^^.rJ  is  obtained. 

In  cylindrical  coordinates  (r,  6,  z),  the  Laplace  equation  V20  =  0  that  is  satisfied 
by  a  potential  function  0  is  defined  by 


r  dr 


f  2^ 
<  dr) 


df\     1  d2<p  d2<p 
'  r2  dd2  +  dz2 


=  0. 


(C.l) 


Replacing  0  by  u^e  wcosn0,  where  K  is  a  real  positive  quantity  and  n  is  an  integer 
number,  results  in  an  ordinary  differential  equation  for  uoo. 


ld_ 
r  dr 


(  du^ 

( 

n2) 

r  — 

+ 

K2 

V  drj 

~r2) 

u  =  0 


(C.2) 


Eq.  (C.2)  is  a  Bessel  equation  of  order  n.  Solution  to  this  equation  may  be  written  in  terms 
of  Bessel  functions  of  first  and  second  kind. 

W(r)  =  AJn(Kr)  +  BYn(Kr)  (C.3) 

Hence,  the  final  solution  to  0  is  obtained. 

0(r,e,z)  =  e'1*  cosnO[AJn(Kr)  +  BYn(Kr)]  (C.4) 

Similarly  a  solution  to  the  Laplace  equation  for  a  potential  function  y/  may  be  written  as 

y/(r.e.z)  =      cos nd\iAJn +  BYnar)]  (C.5) 

where  A  is  a  real  positive  quantity. 
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Green's  theorem  is  expressed  by  following 

JJJ(^V20  -  *V  V>*V  =  JJ(  yr^  -  <P  ^)dS  (C.6) 

where  -—  is  the  directional  derivative  of  <p  in  the  direction  of  the  outer  normal  n  to  the 
an 

surface  element  dS.  The  left  hand  side  of  Eq.  (C.6)  vanishes,  since  both  (j)  and  iff  satisfy 
the  Laplace  equation.  Thus,  Eq.  (C.6 )  is  reduced  to 

where  the  surface  integration  is  being  taken  over  the  surface  bounded  by  a  hollow  cylinder 

a<r  <b  and  the  planes  z  =  0  and  z  =  +°o . 

For  z  =  0,  the  normal  to  the  surface  element  is  expressed  by  n  =  -zeE. 

Consequently  —  at  z  =  0  reduces  to  Kcosnd\AJnacr)  +  BYn(Kr)].  Also  at  z  =  0,  iff  is 
an 


reduced  to  cosnfl 


AJnan+  BYnar)}.  The  substitution  of  y/  and  —  into  the  left  hand  side 


an 

of  Eq.  (C.7),  for  the  flat  end  part  of  the  surface  integral,  yields 


1Kb 


j  j  kcos2  nO[AJnW+  BYn(Kr)][AJnar)  +  BYnar)]rdrdO  (C.8) 

b 

=  7CKj[AJn(Kr)  +  BYn(Kr)^AJnUD  +  BYnar)\rdr. 


0  a 

h 


When  z  =  +°° ,  y/  —  vanishes,  and  nothing  is  contributed  to  the  surface  integral. 
an 

For  the  curved  surface  at  the  inner  radius  r  -  a  the  normal  to  the  surface  element  is 
n  =  -rer  and  consequently       and  y/  may  be  written  as  -Ke'^  cosnd\AJ'n(Ka)  +  BY'n«ai)\ 

and  e"'lzcosn^A7n(Aa)  +  fiyn(^)j,  respectively.  Hence  the  corresponding  part  of  the 
surface  integral  is  defined  by  Eq.  (C.9). 
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jj-Kcos2nde-(K+X)z[Arn(Ka)  +  BY'n(Ka)][AJn(b,)  +  BY^adzdS  (C.9) 

0  0 

=  -f^[AJ^Ka)  +  ^(->JA/.(J-)  +  BYnM] 


For  the  curved  surface  at  the  outer  radius  r  =  b  the  normal  to  the  surface  element  is 
n  =  rer  and  hence,  the  surface  integral  corresponding  to  r  =  b  becomes 


|  J  K  cos2  n6e~{K+X)l[AJ'n(Kb)  +  BY'n(Kb)][Ajn(Xb)  +  BY^bdzdd  (CIO) 

0  0 

= f+^[AJ'n(  rt) + BY'"{  Kb)^J" ( w + **Y"  <;u,)l' 


Working  out  the  other  side  of  Eq.  (C.7)  in  a  similar  way,  and  adding  different  parts  of  the 
surface  integral  together,  yields 

b 

(K-  X)j[AJniKr)  +  BYn^AJnar)  +  BYna,)\rdr  (C-11) 

a 

=  -^—{Ka[AJ'n^BY'n(Ka)\AJn  (Xa)  +  BYn(Aa)\ 

-  Kb[AJ'n(Kb)  +  BY'n(Kb)^AJn<.te)  +  fiyn(AA)] 

-  Xa[AJn(Ka)  +  BYn(Ka)^AJ'n(Xa)+  BY'n(Xa)\ 

+Xb[AJn«b)  +  BYn  (rt)J  AJ'nw  +  BY'n(Xb)^. 

Thus,  if  k  and  A  are  different,  then  we  have 

j[AJn^)  +  BYniKr)j^AJnar)  +  BYnar)]rdr  =  0  (C.12) 

a 

provided  that  <p  and  y/  satisfy  the  homogeneous  Dirichlet  boundary  conditions  for  r  in  the 
interval  a<r<b. 
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[Ajnw+BYn^)]  =  0  [AJHm  +  BYnw]  =  0  (C.13) 

[AJn«a)  +  BYn(Ka)}  =  0  [AJn(Kb)  +  BYn(Kb)]  =  0 

From  conditions  in  Eq.  (C.13)  expressions  for  A,  B,  A  and  B  are  obtained 

A  =  Yn(M>)  B  =  -Jn(Xb)  (C-14) 

A  =  Yn«cb)  B  =  -Jn(Kb) 

where  k  and  A  are  different  roots  of  the  following  transcendental  equation. 

Jn(ifi)YH(ifi)-  J^tpY^tf)  =  0  (C.15) 

Incorporating  the  expressions  for  the  constants  in  the  Eq.  (C.12)  yields 

b 

j[jn(Kr)Yn(Kb)-  Jn(Kb)Yn(Kr)][jnar)Yn(U)-  Jn<,mYnar)]rdr  =  0.  (C.16) 

a 

A  Fourier  Bessel  kernel  based  upon  solution  of  the  Bessel  equation  with  the  homogeneous 
Dirichlet  conditions  in  a  finite  interval  [a,  b]  is  defined  by 

K„iC,.r)  =  [•/„(£/)}; /„(?,<,)}'„ (f.r)]  (C.17) 

where  £,  is  a  root  of  the  transcendental  equation  given  by  Eq.  (C.15).  Introducing  the 
definition  of  ATBc{,.r)  into  Eq.  (C.16)  gives  the  orthogonality  property  of  the  Fourier  Bessel 
kernel. 

b 

jrKni;rr)K^rr)dr  =  0       if  (C.18) 

a 

To  obtain  an  expression  for  the  norm  of  the  Fourier  Bessel  kernel,  we  have  to 
equate  k  and  A  in  the  Eq.  (C.ll).  But,  this  will  induce  both  sides  of  Eq.  (C.ll)  to 
vanish.  However,  we  may  prevent  this  from  happening  by  replacing  K  with  A  +  h  and 
then  by  taking  the  limit  as  h  decreases  indefinitely.  Thus,  Eq.  (C.  1 1)  may  be  written  as 
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b  2  /l  1 

\[Jnar)Yn(Xb)-  Jn(U)Ynar)]  rdr  =  lim(-— — -  (C.19) 

a  1 

{(A  +  h)a[Yn(Xb)J'nia  +  />)«]  -  Jn(Xb)Y'nia  +  h)a]\Y 'n(Xb)J 'n(Xa)  -  J 'n(Xb)Y 

-  (A  +  h)b\Yn(Xb)J'n{(\  +  h)b]  -  Jn(Xb)Y'nia  +  h)b]][Yn(Xb)Jn(Xb)  -  Jn<,xb)Yn(Xb)] 

-  Xa\Yn(mJn[a  +  h)a)  -  Jn(Xb)Ynia  +  h)a]\Yn(Xb)J'n(Xa)  -  Jn(Xb)Y'n(Xa)] 

+  Xb[Yn(Xb)Jn[a  +  hm  -  Jn(Xb)Yn[(x  +      Yn(U>)J'n{»>)  -  J „(^)l^(a*)]}) 

where  A  and  B  have  the  same  expressions  as  A  and  B  respectively,  and  are  replaced  by 
following. 

A  =  A  =  Yn(Xb)  (C.20) 
B  =  B  =  -Jn(Xb) 

Expanding  JHnx  +  h)r]  and  Ynux  +  h)r]  using  Taylor's  theorem  yields  Eq.  (C.21). 

Jn[{X  +  h)r]  =  Jn(Xr)  +  hrJ'n(lr)+0(h2)  (C.21) 
Ynl(X  +  h)r)  =  YnOir)  +  hrY'n(Xr)  +  0{h  ) 

Incorporating  Eq.  (C.21)  into  Eq.  (C.19)  and  taking  the  limit  as  h  goes  to  zero,  results  in 

b  2 

j[JnUr)Yn(Xb)- Jn(Xb)Ynar)]  rdr  (C.22) 

a 

=  ^"{a2  [y„  (Xb)J"(  Xa)  -  Jn  ( AA)F"(Aa)][y„  (Xb)Jn  (la)  -  /„  (Xb)Yn  ( Xa)] 

"I2     „2|"v   ...  i'  .        r  ...v/  .  I2 


+  &2[y„(AA)7n'(^)-  7n(^)yn'(AA)]    -a2[Yn(Xb)J'n(Xa)-  J n(Xb)Y'n(Xa)\ 
"I 


+  ^-\Yn(Xb)Jn(Xa)-  Jn(Xb)Yn(Xa)\Yn(Xb)J'n(Xa)-  J n(Xb)Y'n(Xa)^. 


Since  A  is  a  root  of  the  Eq.  (C.15),  we  may  write 

Jn(xa)Yn(Xb)-Jn(Xb)Yn(Xa)  =  0.  (C.23) 

Incorporating  Eq.  (C.23)  into  Eq.  (C.22)  yields  Eq.  (C.24). 
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0  * 

\[Jn(ir)Yn(U)- Jn(Xb)Ynar)\  rdr  (C.24) 

a 

=  ^b2[Yn(Xb)J'n(Xb)-  Jn(Xb)Y'n(mf  -  a2\Yn{mJ'n(ia)-  J n(Xb)Y'n(Xa)f^ 

Eq.  (C.24)  may  be  simplified  furthermore  by  manipulating  the  terms  in  the  brackets 
of  the  right  hand  side  of  the  equation.  If  y,(D  and  y2(o  are  two  independent  solutions  of 
y"  +  y'f(r)  +  yg(r)  =  0 ,  then  both  satisfy  the  differential  equation. 

yr+tf/('>  +  M('>  =  0  (C.25) 

y?+y'2fir)+y2gw=o  (C.26) 

Eliminating  g(r)  in  Eqs.  (C.25)  and  (C.26)  yields 

(%>f-  yy") + {yiy'\  -  yJi    =  °-  (c.27) 

The  first  parenthesis  in  Eq.  (C.27)  is  the  derivative  of  the  second,  so  that  if  the  latter  is 
denoted  by  z,  Eq.  (C.27)  is  reduced  to  z'  +  zfw  =  0.  Solution  to  this  differential  equation 
may  be  written  in  the  form 

z  =  y2y'l-yyi=Ae~ifl,i*  (C.28) 


where  A  is  an  arbitrary  constant. 

In  the  case  of  the  Bessel  equation,  y"  +  -y'  + 


32  n 

1  ~, 


y  =  0,  we  have  fu)  =  -, 

r 


y,  =  Jnar)  and  y2  =  Ynar).  Thus,  Eq.  (C.28)  becomes 

A[y„ar)7>)-  Y'nar)Jn(Xr)\  =  Ae-lnr  (C.29) 


or 


Ynar)J'nar)  -  Y'nar)Jn(lr)  =  —.  (C.30) 

Ar 
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Weber's  definition  for  the  Bessel  function  of  the  second  kind  K„Ur)  is  given  by  Eq. 


(C.31). 


Y„  (Xr)  = 


COSn7T/n  (Ar)  —  J_n  (Ar) 

sinn7T 


(C.31) 


Substituting  Weber's  definition  of  Ynar),  Eq.  (C.31),  into  Eq.  (C.30)  gives 


.  1  [J'nar)J_nar)~  Jnar)J'_nar)]  = 
smnn  Xr 


(C.32) 


An  expression  for  the  constant  A  is  obtained  by  substituting  the  leading  term  of  the  series 
for  Bessel  function  of  the  first  kind 
Eq.  (C.32). 


°°  (—l)m  (X  r/2)2m+n 
and  Jn(Xr)  =  V  into  the  left  hand  side  of 

m=0      m\T(m  +  n  +  \) 


-1 


sin  n/r 


n  1 


1  -n 


T(n  +  1)  +  1)      r(n  +  1)  IVn  +  1 


J__  A_ 
Xr  Ar 


(C.33) 


71  2 
Since  T(n  +  d  =  nT(n)  and  To  -  n)T<.n)  =  — — — ,  Eq.  (C.33)  reduces  to  A  =  — .  Hence,  Eq. 

71 


smnK 


(C.30)  may  be  written  in  the  form 


Y„(lr)J'n(Xr)-  YjXr)Jn(Xr)  = 


n v    '  n  * 


-2 
n:Xr 


(C.34) 


Evaluating  Eq.  (C.34)  at  r  =  b  results  in 


YniXb)J'n(Xb) -  Y'n(ib)Jn(m  = 


-2 
nXb 


(C.35) 


Evaluating  Eq.  (C.34)  at  r  =  a  and  using  Eq.  (C.23)  yields  Eq.  (C.36). 

-2   J  (lb) 


Yn(Xb)J'naa)-Y'„(Xa)J„{Xb):= 


7tXa  J„(Xa) 


(C.36) 


We  may  rewrite  Eqs.  (C.35)  and  (C.36)  in  terms  of  Kn((rr). 
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Br 


2  Jn<,(,h) 
Jia  /„(f,.<i) 


dK 


Sail 


dr 


r=b 


_2_ 
Kb 


(C.37) 


Incorporating  the  expressions  from  Eqs.  (C.35)  and  (C.36)  into  Eq.  (C.24)  gives 

L  9 

2 


j[jn(Xr)Yn(u,)-Jn(U)Ynar)]  rdr  =  -r^ 


(C.38) 


By  rewriting  the  left  hand  side  of  Eq.  (C.38)  as  ||/Cn(f  ,r)|,  the  expression  for  the  norm  of 
the  Fourier  Bessel  kernel  is  obtained  from  Eq.  (C.39). 


v      II2  2 


K  £ 


1- 


(C.39) 


APPENDIX  D 
STRESSES  IN  TERM  OF  STRESS  FUNCTIONS 


In  this  appendix  stresses  are  derived  from  the  axial  and  radial  stress  functions. 
The  displacement  vector  is  defined  in  terms  of  stress  functions  by  Eq.  (D.l). 

u  =  -aV2(fr(r,z)er+  (Pz(r,2)ez)  +  («  +  i3)V[V.('Pr(r.2)er+  fz(r,z)ej  (D.l) 


Rewriting  Eq.  (D.l)  in  the  component  form  yields 


u,  =  -a 


r  J 


uz  =  -aV2f2 +(<*  +  /?)  . 

dz 


dr 


(D.l) 


where  q>  is  given  by  Eq.  (D.3). 


9  =  -— (r«Pr)  +  — -i 
r  dr  dz 


(D.3) 


From  the  constitutive  equation  and  strain-displacement  relationship,  the  stress  tensor  may 
be  written  in  the  form  of  displacement  vector. 


a  =  A  (V.u)I  +  ^[Vu  +  VTu] 


(D.4) 


By  using  Eq.  (D.4)  stresses  are  obtained  in  terms  of  displacements. 


a  =  X 


Id,  .  du 
r  dr  dz 


+  2/1 


du 


(D.5) 
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ae  =  X 


Id  .  .  du 
r  dr  dz 


r 


(D.6) 


Id.  .  du7 
r  dr  dz 


_  du? 

^Tz 


(D.7) 


'  dur     du,  \ 
— -+ — L 

V  dz  dr 


(D.8) 


Introducing  Eq.  (D.2)  into  Eqs.  (D.5)  through  (D.8)  gives, 


.    1  d 
or  =  -Aa — — 
r  dr 


V— U 

r 


aA  d{  dq>\ 
r  dry   dr  J 


(D.9) 


dz  dz\dz  J 


-2/ia-l  V2- 


+  2//(a  +  /?) 


dr 


v,  dr  j 


.     1  5 
cre=-Aa-— 
r  or 


V2- 


r2y 


•p. 


+  A(a  +  /J)- 


\d(  dq>^ 


dr 


V   dr  j 


(D.10) 


-A«|(V2»Pz)  +  A(a^)|-f^ 
e>Z  dz\dz  ) 


-2/ia 


.2 


V 


r  y 


+  2//(a+/3)- 


J      1  * 

o\  =  -A  a — — 
r  <9r 


v2-4i 


+  A(a  +  /?)-— 
r  dr 


I  d  (  d(p 


dr 


V    or  j 


(D.ll) 


dz  dz\dz ) 


dz 


+  2n(a  +  p)— 


(        1  A 
r  ) 


(D.12) 


The  terms  — — 
r  or 


y2  — rVr  and  -T-(v2,pJ  may  be  written  as  V: 


(r«F,) 


and 


—  f 


,  respectively.  Hence,  substitution  of  these  terms  into  Eqs.  (D.9)-(D.  1 1)  leads 


to 


a.  =  -A  aV' 


1  d  /  f<?*o 


-2\ia 


{  r2) 


dr 


dz 


+  X(a  +  p) 


r  dr 

d_(d(p^ 
dr 


\_d_f  dtp\  d2(p 
V  dry  dz2 


(D.13) 


\dr  j 


oB  =  -A  aV' 


+  A(a  +  /3)-— 
r  e>r 


1  <9  (  dip} 


V 


-AaV: 


dz  r 


v  41 

T  2 

V      r  ; 


•Fr  +  2M«  +  i3)-f<9<i£) 
r 


•P. 


(D.14) 


o,  =  -A  aV3 


+  A(a  +  )8) 


1  d 


f  dcp^ 


+  X(a  +  p)^-2^iaV2(^- 

dz  \dz 


r  dr 


V  Br  j 


+  2/i(a  +  P) 


-AaV2 
d2cp 


( d  ^ 
^dz 


(D.15) 


J 


dz2' 


1  (9 


<9f,  . 


1  ^ 


'  dq>\    d2q>  . 


but  — — (rfr)  +  — is  p  and 


V  dr)  dz2 


is  the  Laplacian  of  (p.   Thus,  Eq. 


(D.  13)-(D.  15)  are  reduced  to 


ar  =  A/3V>-2//a— 
dr 


Wr  +  2n(ct  +  p) 


d2cp 
dr2 


(D.16) 


a,  =  A  pv2cp  -  2»al-h2  - ±  V  +  2/*(a + 0)±f  ^ 
r  V       r  V  r\dr  j 


(DAI) 
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a  =A^V>-2//aV: 


( d  \ 


+2n(a+p) 


d2cp 
dz2' 


(D.18) 


By  replacing  a  with  -— ,  B  with  — - —  and  X  with    2^V  ,  Eqs.  (D.16)-(D.18) 

A  +  2^  l-2v 


transform  to 


a  1  v2  -  — 1 


V       r  y 


1-  v 


)2  A 


vV2 


dr2 


<P 


(D.19) 


0o  = 


2 

\        r  ) 


1      f    V72  1 

«P  +   vv2--— 

r  or  ) 


1-  v 


(D.20) 


vV2- 


(p. 


(D.21) 


The  same  substitution  into  Eq.  (D.12)  (shear  stress)  results  in, 


d2q> 


!P  — 

1  -  V  drdz 


(D.22) 


Eqs.  (D.19)-(D.22)  define  the  stress  components  in  terms  of  stress  functions. 


APPENDIX  E 

STRESSES  IN  TERMS  OF  FOURIER  BESSEL  KERNELS 


In  this  appendix  stresses  are  formulated  in  terms  of  Fourier  Bessel  kernels  and 

trigonometric  functions. 

Stress  functions  Jfyr,*)  and  P^r,*)  are  defined  by  I    I  Xr(r'.z')Grir,r',z,i')drdz 

and        X<r',z')GAr,r\z,z')dr'dz'  where  Gr(r,r-,z,z')  and  G,(r,r-,z,z')  are  radial  and  axial 

J-LJa     z  z 

Green's  functions,  respectively,  given  by  Eqs.  (E.l)  and  (E.2). 


Grir,r-,z,z-)  =  ~2,i  —  —2  ~2  rKx(^.r)Kx 


n/rz  «;rz 
■cos  cos 


(E.l) 


Gz(r,r',z,z')  =  -  —  ^^ 
L    (=1  n=l 


Co,-^Q(Co,-a) 


(L2Co2,+n2^2)2/o2(^)-yo2(V) 


„  nnz  .  nnz 

r  ^0(fw.r-)^0(fw,r)Sin— sin— 


(E.2) 


From  Eqs.  (E.l)  and  (E.2)  and  by  introducing  new  functions  Tr^u,n)  and  rz(£0.,n),  the 
stress  functions  may  be  written  as 


oo  oo 


^('.d  =  yy  rr(fu ,  n)Kx  (Cu .  0  COS 


n/rz 


,=l  n=0 


(E.3) 


/=1  n=l 


(E.4) 


where  rr(fli>n)  and  rz(?0i.„)  are  given  by  Eqs.  (E.5)  and  (E.6). 
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!%((„..„)  = 


rr(f„,o)  =  - 


K2  L 
7T2         Z.2^)  1      fL  f*  /V         I-  j  /  i  / 


(E.5) 


Co^o(Co,") 


(L2Co2,+nV)27o(?„f<.)-702(V, 


(E.6) 


The  components  of  stress  in  terms  of  stress  functions  are  given  by  Eqs.  (E.7)-(E.10), 
where  in  these  equations  (p  is  given  by  Eq.  (E.l  1). 


a  =2- 


dr 


V        r  J 


1 


1- V 


vV2- 


(E.7) 


v     r  y 


r  1-v 


r  dr) 


(E.8) 


1-v 


)2  A 


vV2-- 


A2 


9> 


(E.9) 


V       r  y 


r    1-v  <9r<9z 


f  - 


(E.10) 


<P  =  -T-(A-^r)  +  -V- 


(E.ll) 
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The  substitution  of  Eqs.  (E.3)  and  (E.4)  into  Eqs.  (E.7)-(E.10)  gives 


oo  oo 


i=l  n=0 


(v-4) 

V        r  ) 


J^(f,.,r)COS 


1      ^  M  n=0 


'       2  ^Mia 


oo  oo 


vV- 


I  —  V  i=\  „=o 


dr 


vV2- 


J 

d2\d 


dr 


rK.(C,,r)COS 


nnz 


dr2 


dz 


tf0(f„„--)sin 


nnz 


(E.12) 


°,»=2XXr'(f»',)' 

1=1  n=0 

1  OO  oo 

+rLHr 


"V-41 

2 


v 


<£„.«) 


#.(£.,,-)  COS 


n?rz 


V72      1  d\ 

VV  — 

r  dr  J 


\_d_ 
r  dr 


rAT,  (£,.,;•)  cos 


AlTTZ 


1  <9  ^  <9 


vV2--^- 
r  dr  J 


dz 


#0(C0,,r)Sin 


/zttz 


(E.13) 


/=1  n=0 


Aocc.Dsin 


«;rz 


1     v  i=\  «=o 

1  oo  oo 


vV2 


L 

<9z2  J  r  <9r 


rA^.ocos 


n^z 


vV2- 


dr2 


dz 


'  .  nKz~\ 


(E.14) 


1  =  1  n=0 


„  nnz 
K0((v.r)sm— 


1  =  1  n=0 


rt7TZ 


i  <?2  n 

-7— ;2,2,rr««.-> 

i        V  I=1  n=0 


<9r<9z 


dr  a. 


rK^c  u,r)  cos 


n^z 


1       V  ,=1  „=0 


<9r<9z  cfe 


^0(fo,-r)sin 


rc/rz 


(E.15) 
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The  Fourier  Bessel  kernel  ATj(fu,r)  is  obtained  from  Bessel  functions  of  order  one. 
Hence,  it  satisfies  Bessel's  equation  which  translates  to 


r  dr\  dr )  r2 


Kxa;h.r)  =  -QiKx^rr). 


(E.16) 


n 


Also  from  basic  identities  for  Bessel  functions,  J^cur)  =  Ju  ,(f1(r)  ^(f,/)  and 

Si/ 


fyf,/)  =  y„-><f,/>  -  we  obtain 


(9  1 

—  K,  (fB.r)  =  &j[  J0«^>]j (f„W  -  /o<C» yi  <V>]  ^1  <&u ■ r)- 

dr  r 


(E.17) 


Likewise  for  K0^orr)  we  may  write 


ld(rd\ 


r  dr\  dr 


(E.18) 


dr 


(E.19) 


Using  results  of  Eqs.  (E.16)-(E.19)  and  introducing  new  kernels  Kwau.r)  and 
K]0(torr),  Eqs.  (E.12)-(E.15)  are  reduced  to 


00  00 


(£„■.") 


i=l  n=0 

1      '  i=l  n=0 
»  1       °°  ~ 

1- vtr^ 


r  ,2  »vy.  „  iv 


cos- 


nnz 


2  _2  \ 


n  7T 


v 


-Cu  7T  Ci«*oi<Cw'>cos 


n7rz 


I 


V  r  J 


cos- 


<Cw«> 


1  ~  V  i=l  n=0 


L  ) 

Tj 


-c2  -■ 


,2—2  \ 


#0«;0i,r)cos 


2     n  7T 

f  1 

V  r 


nnz 


cos- 


(E.20) 
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oo      oo  |  I 

^=2xirr(c,,n,-  -a- 


1=1  n=0 


if,  (C...r) COS 


n7K 


oo  oo 


1  -  V  j=,  fl=0  V  L 


„2_2  V 


CA«..')COB 


I  oo       oo  I 

1      "  (=1  n=0 


r)COS 


1  oo       oo  Yl'JL 

1  -  V  ,=|  n=0  l  V 


oo  oo 


~I3 


n7Tz 

2^-2  > 


r)COS 


n7Tz 


nn  1 


1      V  ,=i  n=o  ^ 


«7EZ 


(E.21) 


oo  oo 


nnf  n2n2^ 

°z  =  2LLTz^-r  -Co,-  — 7T 

■*■>]  „=o  *-  V  L  J 


K0((orr)COS 


nnz 


oo  oo 


n  n 


i    v  ,=1  „=0 


C„^0,(fu,r)COS 


nTEz 


DCOS 


nnz 


i'=l  n=0 


1 


Co,^O^o,-r)--^10^o,-r) 

r  j 


cos- 


n^z 


(E.22) 


oo  oo 


n2n2^ 


V  L 


1=1  n=0 


«  „  nnz 
-Co/^io  ««■'■)  Sin-T" 


+I£rr(fll,.)H 


i=l  n=0 


„2_2  \ 


nnz 

tfX.osin— - 


i      v  ,=1  n=0 
1 

"irrXSr 


UK 
~L 
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where  K0l {(u , o  and  Ki0(Cv,r)  are  given  by  Eq.  (E.24)  and  (E.25). 
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Note  that  to  derive  Eqs.  (E-20)-(E.23),  the  following  identities  were  used. 
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Eqs.  (E-20)-(E.23)  are  simplified  further  by  combining  the  summations  and  canceling  the 
common  terms. 
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Eqs.  (E.26)-(E.29)  describe  the  components  of  stress  in  terms  of  Fourier  Bessel  kernels 
and  trigonometric  functions. 


APPENDIX  F 

EQUATIONS  FOR  THE  ARBITRARY  CONSTANTS  OF  THE  COMPLEMENTARY 

SOLUTION 

In  this  appendix  radial  boundary  conditions  for  £o-,z)  are  applied  and  four  linear 
equations  for  the  arbitrary  constants  of  £(r,o  are  derived 

AAA  /V 

Expression  for  %u,z)  is  given  by  Eq.  (F-l),  where  A,-,  \,  Bn,  Cn  and  Dn  are 
arbitrary  constants. 
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Radial  boundary  conditions  for  §(r,*)  are  given  by  Eqs.  (F.2)-(F.5). 
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where        through  p4(«)  are  known  functions. 

From  Eq.  (F.l)  and  the  definitions  for  the  Fourier  Bessel  kernels,  we  may  write  the 
following. 
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The  substitution  of  Eqs.  (F.13)  and  (F.10)  into  radial  stress  boundary  condition  at 
r  =  a ,  yields 
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Using  the  orthogonality  property  of  cosine  and  introducing  new  functions  A,  (£„,«),  X, ,(,.), 
X12(«>,  X13(n)  and  X,4(«),  Eq.  (F.14)  reduces  to  Eq.  (F.15). 
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The  definitions  for  X,,(n),  X,2(n),  X13(«),  X14(«)  and  A,{fof.«)  are  given  by  Eqs.  (F.16) 
through  (F.20),  respectively. 
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Where  in  Eq.  (F.20),  Q^.n)  and  Q2(?o.-">  816  given  by  Eqs.  (F.21)  and  (F.22). 
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The  substitution  of  Eqs.  (F.6)  and  (F.12)  into  the  radial  stress  boundary  condition 
at  r  =  b ,  yields 
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where  X21</o,  X22<n),  X23<n),  X24u)  and  A2(col,«)  are  given  by  Eqs.  (F.24)  through  (F.28), 
respectively. 
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The  substitution  of  Eqs.  (F.7)  and  (F.9)  into  shear  stress  boundary  condition  at 
r  =  a ,  yields 
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where  X,,(«),  X„  and  A3(foj.n)  are  expressed  by  Eq.  (F.30)-(F.34). 
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The  definitions  for  £23(f0i,«)  and  £24  (?„,..«)  are  given  by  Eqs.  (F.35)  and  (F.36). 
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The  substitution  of  Eqs.  (F.8)  and  (F.ll)  into  the  shear  stress  boundary  condition 
at  r  =  b ,  yields 
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where  X4,(«),  X42(«),  X43(«),  X^u)  and  A4(foj,«)  are  defined  by  Eqs.  (F.38)-(F.42). 
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Eqs.  (F.15),  (F.23),  (F.29)  and  (F.37)  represent  four  linear  equations  for  the  five  arbitrary 
constants  A/,  An,  Bn,  Cn  and  Dn . 
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